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ABSTRACT 


The  spatial  distribution  of  power  in  a  nuclear 
reactor  tends  to  follow  the  neutron  distribution  because  the 
energy  resulting  from  nuclear  fission  is  released  largely 
at  the  point  where  neutrons  induce  fission.  Leakage  in  an 
unreflected  reactor  causes  both  distributions  to  be  higher 
near  the  reactor  center.  The  power  at  the  center  is  limited 
by  the  maximum  temperature  which  can  be  tolerated  by  core  ma- 
terials with  the  result  that  outer  regions  operate  at  less 
than  ultimate  capacity.   Ideally  all  regions  should  be  oper- 
ated at  the  maximum  tolerable  temperature,  which  would  require 
uniform  spatial  distribution  of  power.   Thus  an  important 
improvement  in  a  reactor* s  power  output  can  be  achieved  by 
making  the  power  distribution  more  nearly  uniform. 

One  simple,  economical  method  of  achieving  uniform 
power  consists  of  compensating  for  the  decrease  in  power  near 
the  edge  of  the  core  by  means  of  a  reflector.   If  the  return 
of  neutrons  from  the  reflector  to  the  core  can  be  adjusted  to 
compensate  for  the  leakage  factor,  the  power  distribution  can 
be  markedly  improved.   This  method,  applied  to  a  reflected 
spherical  homogeneous  reactor,  is  studied  in  this  paper. 

The  distribution  of  neutrons  in  a  reactor  is 
specified  reasonably  well  by  the  diffusion  type  equations  to 
be  discussed  in  Chapters  I  and  II.  Multi-group  methods  of 
solution  are  used  because  of  the  importance  of  the  dis- 
tribution of  fast  neutrons.   Assuming  two  groups  of  neutrons 
was  considered  inadequate  for  the  intermediate  type  of  re- 
reactors  which  are  studied.   Multi-group  calculations  are  so 
time-consuming  that  the  use  of  computing  machines  is  essential, 
Chapters  I  and  II  develop  techniques  for  use  with  M.  I.  T.'i 
digital  computer  Whirlwind  I. 


ABSTRACT 
(continued) 

To  simplify  the  problem,  a  basic  reactor  is 
selected  with  highly  enriched  fuel,  sodium  coolant,  beryl- 
lium moderator  and  reflector.  The  reactor  is  assumed  to  be 
homogeneous  and  in  the  steady  state.  The  effect  of  reactor 
poisoning  is  not  considered.  Spherical  geometry  is  assumed. 

It  was  found  that  by  suitable  adjustment  of  the 
material  proportions  in  both  core  and  reflector  the  average 
power  could  be  raised  to  within  16%  of  the  maximum  power. 
Results  are  presented  and  discussed  in  Chapter  III. 
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IPO.-.;  POWER  DISTRIBUTION  IN 
FLECTED  KO     "OoC  REACTOR 

Neutron  leakage  from  an  unreflected  reactor  results 
in  a  lower  concentration  of  neutrons  toward  the  reactor  edge. 
By  virtue  of  the  fact  that  the  power  resulting  from  nuclear 
fission  induced  by  neutrons  is  largely  liberated  at  the 
point  where  the  fission  occurs,  the  spatial  distribution  of 
power  also  is  depressed  near  the  reactor  edge.  This  imposes 
a  limitation  on  a  reactor's  power  output.  The  maximum  power 
density  at  any  point  is  limited  by  temperature  considerations. 
Outer  regions,  with  lower  power  density,  operate  at  tem- 
peratures lower  than  their  ultimate  capacity.  Elementary 
considerations  indicate  that  the  average  pov/er  in  an 
unreflected  spherical  reactor  is  only  one  third  of  the 
maximum  power  produced  at  the  reactor  center.  Thus,  if  a 
uniform  power  distribution  could  be  achieved,  an  important 
increase  in  reactor  output  could  be  realized. 

Spatial  variations  of  the  proportions  of  the  mate- 
rials in  the  core,  together  with  use  of  a  reflector,  is  one 
method  of  improving  pov/er  distribution.  Another  is  the  use 
of  a  poison  to  depress  the  neutron  concentration  near  the 
reactor  center.  These  methods  have  the  disadvantages  of 
increased  complexity  and  cost.  A  simple,  economical  alter- 
native, the  use  of  neutrons  reflected  back  into  the  core 
from  the  reflector  to  compensate  for  the  power  depression  in 
the  outer  core,  is  investigated  here. 


The  diffusion  equations  to  be  discussed  in 
Chapters  I  and  II  specify  reasonably  well  the  distribution 
of  neutrons  In  a  reactor.   Since  neutron  energy  varies 
widely  and  since  the  coefficients  of  the  terms  in  the  dif- 
fusion equation  depend  on  neutron  energy,  some  simplifi- 
cation is  necessary.   One  method,  applied  to  thermal  re- 
actors, assumes  that  neutrons  are  slowed  to  thermal  energies 
immediately  after  birth.   The  coefficients  of  the  terms  in 
the  diffusion  equation  are  then  determined  at  thermal 
energies,  and  a  simple  analytical  solution  is  possible. 

However,  in  this  problem,  the  number  of  neutrons 
with  energies  greater  than  thermal  must  be  increased  in  order 
to  supply  the  reflector  with  fart  neutrons  for  moderation  end 
reflection  back  into  the  core  to  compensate  for  leakage. 
Therefore,  the  distribution  of  higher  energy  neutrons  must 
be  known,  and  one-group  calculations  are  inadequate. 

Considering  two  groups  of  neutrons — fast  and 
slow — would  give  some  representation  of  fast  neutrons  whioh 
leak  from  the  core  to  the  reflector,  where  they  may  be 
moderated  and  returned  to  the  core  as  slow  neutrons.   Hew- 
ever,  better  solutions  are  possible  if  more  than  two  groups 
are  considered. 

Multi-group  calculations  require  solutions  of  dif- 
fusion equations  for  each  energy  group  of  neutrons  con- 
sidered and  are  very  laborious  and  time-con3uming.   The  use 
of  a  high  speed  digital  computer  serves  to  minimize  both  of 
these  problems.   In  Chapters  I  and  II,  diffusion  equations  foi 


four  and  six  energy  groups  respectively  are  broken  down  Into 
a  set  of  simple  calculations  which  a  computer  can  perform  to 
solve  the  equations.   The  actual  programmed  set  of  in- 
structions is  used  with  M.  I.  T.fs  digital  computer, 
Whirlwind  I.   The  basic  method  is  applicable  to  other 
digital  computers  by  a  suitable  transformation  of  the 
specific  set  of  instructions.   Whirlwind  I  requires  less 
than  a  minute  to  solve  the  equations,  but  a  period  of  from 
two  to  three  minutes  is  required  for  the  computer  to  trans- 
form the  input  of  decimal  numbers  and  interpreted  instruc- 
tions to  its  own  binary  system.   Thin  conversion  time  can  be 
reduced  by  rewriting  the  program  in  the  basic  Whirlwind  code, 
The  Whirlwind  I  computer  will  perform  this  function  auto- 
matically and  will  cut  a  new  tape  containing  the  converted 
program. 

In  this  paper  a  basic  reactor  composition  and  con- 
figuration are  chosen,  and  the  proportions  of  the  various  ma- 
terials of  the  core  and  reflector  varied  in  order  to  de- 
termine their  effect  on  the  power  distribution.   The  funda- 
mental reactor  type  studied  was  one  considered  suitable  for 
high  output  mobile  reactor  application. 

The  requirements  of  this  type  of  reactor  led  to 
the  use  of  core  materials  capable  of  withstanding  hi,gh 
temperatures.   A  liquid  metal  coolant  was  chosen  in  order  to 
avoid  having  to  pressurize  the  core.   Of  the  liquid  metals 
available,  sodium  was  chosen  because  of  its  superior  heat 

ill 


removal  characteristics  and  because  of  its  compatibility  with 
other  materials  at  high  temperature  with  respect  to  corrosion. 

Once  sodium  was  chosen  as  the  coolant ,  stainless 
steel  became  the  best  material  for  cladding  in  view  of  its 
high  temoerature  strength  and  resistance  to  corrosion  by- 
sodium. 

A  dispersal  fuel  element  of  uranium  dioxide  in  a 
matrix  of  stainless  steel  was  chosen  to  meet  high  temperature 
requirements  and  to  minimize  radiation  damage.   Calculations 
for  a  fuel  element  of  this  type  have  been  made. 

Beryllium  was  chosen  as  a  moderator  and  coolant 
because  of  its  high  moderating  ratio  and  good  thermal  proper- 
ties. 

The  average  operating  temperature  of  the  reactor 

o 

was  taken  as  1000  F,  which  corresponds  to  a  thermal  neutron 

energy  of  0.07  ev. 

For  mathematical  simplification,  spherical  geometry 
was  assumed.   The  reflector  thickness  was  uniform.   Both  core 
and  reflector  were  assumed  to  be  homogeneous  and  in  the 
steady  state.   The  effect  of  reactor  poisons  was  not  con- 
sidered.  Thus  an  actual  reactor  would  require  additional  re- 
activity to  overcome  the  effects  of  fission  product  poisons. 


"Proceedings  of  the  International  Conference  on  the  Peaceful 
Uses  of  Atomic  Energy,"  Volume  9>  page  561,  C.  E.  Weber  and 
H.  H.  Hlrsch. 
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CHAPTSH  I 
Method  of  Finite  Differences 
L.  N.  Hoover 

This  chapter  will  develop  a  program  for  solving 
multigroup  neutron  diffusion  problems  in  spherical  geometry 
on  the  Whirlwind  I  (WWI)  digital  computer  using  a  finite 
difference  approximation  to  the  age  diffusion  equation.   The 
solution  will  give  the  power  distribution,  the  neutron  flux 
distribution  in  eaoh  of  several  energy  groups,  and  the 
critical  radius  for  a  reflected  spherical  nuclear  reactor 
having  a  homogeneous  core  and  a  homogeneous  reflector.   The 
program  will  be  developed  for  neutrons  in  four  energy  groups. 
By  suitable  modification  of  the  ba3ic  program,  the  number  of 
energy  groups  and  the  number  of  regions  considered  can  be 
increased.  Four  energy  groups  represent  the  maximum  number 
of  groups  which  this  program  will  conveniently  handle  using 
only  the  internal  storage  of  WWI.  To  increase  the  number  of 
energy  groups  considered  will  require  using  auxiliary 
storage,  whloh  has  the  disadvantage  of  increasing  the 
machine  time  required  for  a  solution.   The  four-group 
program  requires  approximately  three  minutes  of  computer 
time  for  a  complete  solution.  For  a  parameter  study,  the  use 
of  a  larger  number  of  energy  groups  was  not  considered  to  be 
advantageous,  since  very  little  improvement  in  results  would 
be  achieved  by  the  use  of  more  groups. 


I.   ^ge  Diffusion  Equation 

The  age  diffusion  equation  provides  a  good  de- 
scription of  the  neutron  behavior  of  a  large  class  of 
nonthermal  reactors.  The  following  development  is  concerned 
with  reactor  assemblies  that  are  large  compared  with  the 
neutron  mean  free  path,  and  in  which  the  neutron  energy  is 
reduced  primarily  by  elastic  collisions.  The  method  accounts 
for  both  fact  and  thermal  fission,  for  the  fission  spectrum, 
and  for  inelastic  scattering. 

The  basic  equation  used  is  the  age-cliff usion,  which 

gives  good  results  for  the  case  of  weak  capture  and  a 

2 
cross  section  that  varies  slowly  with  energy, 

-DV20(u)  +  f.  (u)0(u)  =  "  ^U}   *  S(£,U)     (1.1) 

The  term  -  — e-* — '—   represents  the  excess  neutrons  entering  the 
a    u 

lethargy  interval  over  those  leaving,  assuming  that  the  slow- 
ing down  density,  q,  can  be  treated  as  a  continuous  function 
of  the  lethargy,  u.  This  excess  number  of  neutrons  is 
balanced  in  the  steady  state  by  the  leakage  of  neutrons  out  of 
the  region  on  interect,  which  is  -DV  0(u),  where  0(u)  is  the 
flux  per  unit  lethargy  interval,  and  by  the  neutron  absorption, 

2-Q(u)  0(u).  The  term  S(r,u)  represents  neutrons  entering 
a 

lethargy  u  at  a  space  point  r  as  a  result  of  production  by 
fissions,  plus  inelastic  scattering  from  a  higher  energy. 

There  are  advantages  to  expressing  Equation  (1.1) 
in  terms  of  the  slowing  down  density,  q,  rather  than  the 
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flux,  0»  The  relationship  between  the  flux  per  unit  lethargy 

interval  and  the  slowing  down  density  in  the  case  of  slowly 

2 

varying  capture  cross  section  and  weak  absorption  is 

0(u)  Z0V  (1'2) 

For  energies  below  the  Bamoe  energy,  the  same  number  of 
neutrons  must  slow  down  past  each  energy,  except  for  the  ef- 
fects of  leakage  and  absorption.   Hence,  for  reactors  with 
small  absorption  in  the  slowing  down  range,  the  slowing  down 
density  tends  to  be  a  less  sensitive  function  of  energy  than 
the  flux,  0,   which  will  vary  more  rapidly  than  q  as  2.3(u) 
varies* 

In  terms  of  q,  Equation  (l.l)  becomes 


-D\72 


fXs(u>  J  j2ju>         d   u      -  (1#; 


In  Equation  (1.3),  u  is  a  continuous  variable* 
For  multigroup  calculations,  Equation  (1*3)  is  Integrated 
over  the  lethargy  range  from  U]_  to  Mot   corresponding  to  the 
width  of  the  i'th  energy  group,   Z!  g(u)  is  considered  as 
Independent  of  position  and  brought  in  front  of  the 

Laplacian  operator, V  •   The  age- diffusion  equation  for 

4 

lethargy  group  i  then  becomes 

Tr  (  D  __2J        T7  (j[a  q)  B  OUT 

i^s    ji  avg     i^)*-8  '1  avg 

+  qiIN  +   js^ujdu,  ix.k) 
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u2 


v/here  \)*    =  /    du.  The  averages  are  with  respect  to  u  in 

1      IN 
the  group.   q.    is  the  degradation  of  neutrons  into  the 

i'tii  group  from  the  group  of  next  higher  energy,  and  q^     is 
the  degradation  out  of  the  i  'th  group.  The  approximation  of 
replacing  the  averages  of  products  appearing  in  Equation(1.4) 
by  the  product  of  averages  yields 

->   '  1  avg  -5   s  i  avg       u 

+  ^i°UT  =  ^i™   +  /  S(r,u)du         (1.5) 

/Ui 

This  procedure  is  justified  if  either  q  or  the  cross  sections 
do  not  change  radically  over  the  group. 

To  solve  Equation  (1.5)  when  q.  '  and  S  are  known, 

it  is  necessary  to  postulate  an  additional  relationship 

between  (q)j_  avc;,  q.%      ,    and  q$_        •  The  simplest  method  would 

be  to  assume  a  linear  relationship  for  q  across  the  lethargy 

group  and  to.  take  an  arithmetic  average.  However,  with 

only  four  lethargy  groups,  the  large  range  of  lethargy  per 

group  tends  to  make  the  linear  relationship  inaccurate.  A 

suitable  method  for  use  with  a  small  number  of  energy  groups 

is  to  use  the  relationship 

(o^i  avg  =  fcU  (1*6) 

where  f  is  a  number  chosen  by  trial  and  error  for  each 

group  and  region  of  the  assembly  so  that  f  times  the  spatial 

OUT 
average  of  qj_    gives  a  correct  value  for  the  spatial 

average  of  (q)j_  aVrr»  The  solution  is  not  sensitive  to  the 

choices  of  f,  and  it  is  usually  possible,  on  the  basis  of 
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experience,  to  make  adequate  choices  the  first  time.  From 
the  problem  solution,  the  values  of  (q)^  aVe  are  plotted 
versus  lethargy  as  in  Figure  12,  A  smooth  curve  is  faired 
through  these  values  of  (q)i  av~  to  represent  the  actual 
process  of  continuous  slowing  down.  From  this  curve, 
refined  values  of  f  can  be  determined  and  used  to  improve  the 
solution,  if  the  initial  estimate  was  greatly  In  error. 
II.  Average  Nuclear  Properties  of  Material  in  a  Lethargy 
Group 

Equation  (1.5)  expresses  the  age  diffusion  equation 
for  lethargy  group  1  in  terms  of  material  properties  averaged 
over  the  i'th  group.  Within  each  group,  the  neutrons  are 
assumed  to  diffuse  without  energy  loss  until  they  have  suf- 
fered the  average  number  of  collisions  which  would  be 
required  to  decrease  their  energy  to  that  of  the  next  lower 
group.  At  this  point  it  is  supposed  that  the  neutrons  are 
suddenly  transferred  to  the  latter  group.  This  process  is 
supposed  to  continue  while  the  neutron  energy  is  degraded 
from  the  group  of  birth  energy  to  that  of  lowest  (thermal) 
energy. 

Since  the  energy  spectrum  of  the  neutrons  in  a  re- 
actor covers  a  continuous  range  from  thermal  energies  to 
about  10  Mev,  the  group-diffusion  concept  is  evidently  an 
approximation.   This  is  partially  overcome  by  assigning  ap- 
propriate average  values  for  each  group  to  the  various 
properties,  such  as  the  cross  section  and  the  diffusion 


-6- 

coefficient  of  the  core  and  the  reflector. 

a.  Diffusion  Coefficient.  The  average  diffusion 
coefficient,  Dj_ ,  for  energy  group  i  is  derived  from  a  con- 
sideration of  the  fust  neutron  current  density.  The 

3 
definition  of  Dj_  is 

J^1D(E)^70(r,E)dE 

Di    =    T^T (2-1} 

/    V  0(r,E)dE 

It  can  "be  seen  from  this  equation  that  D.  has  the  same  value 
at  different  spatial  points  within  the  region  being  con- 
sidered only  if  0(r,E)  is  separable  into  a  function  of  r 
multiplied  "by  a  function  of  E,  i.e.,  if  the  variation  in  the 
neutron  flux  spectrum  from  point  to  point  in  the  medium  can 
be  neglected.  This  assumption  is  usually  made  for  simplicity. 

Hence,  with  D(E)  replaced  by  1 ,  the  expression  for  D* 

3ftr(E) 
becomes        /"En     -, 


E2   '3ltr(E)  0(E) dE 


;E)dE 


/  Eo 


For  the  asymptotic  case  where  E«oCEQ,  which  is  assumed  to 
apply  for  all  energy  groups  considered  here  using  Be  as 
the  moderator,  and  for  the  case  of  weak  absorption 

0(E)  «   *W  ,   or0(u)  ="c^k 
In  lethargy  notation  we  find  that  Equation  (2.2)  becomes 


-7- 

...   q(y>l-dwl  . 
*qU)  fo 


u1  (     q     )  (       i     ) 

Ui<fz7>lavg  Wlavgqi^ 

The  last  step  Involves  the  approximation  of  replacing  the 
averages  of  products  by  the  product  of  averages.  This  is 
Justified  if  either  q  or  the  cross  sections  do  not  change 
radically  over  the  group.  Hence,  this  is  a  criterion  which 
must  be  considered  in  choosing  the  energy  subdivisions  of 
the  group.  Thus  the  average  diffusion  coefficient  for 
group  i  is  given  by 

(    1   ) 
Di  z    <3jggrtr)1  ftYg  (2oJ 

(  e?  ) 

j  ^s  i  avg 

b.   Total  Flux  Per  Energy  Group*  The  total  flux 
in  energy  group  i,  NVj_ ,  is  defined  as 


NV±  -   )0(u)du=    q(u)du 


u: 


Use  of  the  same  approximations  as  for  the  diffusion  coef- 
ficient  gives  the  result,  where  17^  -  J    du 

NV1  =  Ull7Wi  avg  U)l  «B  ^'^ 


or 
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NV1 


^h   avg  =    <   2  j  (2.^a) 

"j^B^l   avg 

c •   Average  Absorption  Crp s 9  Section  Per  Group. 

The  average  absorption  cross  section  for  energy  group  .1  Is 

defined  as 

"a 

r   _  -4£a(u)0(u)du  -  u,  $z^*\ 

When  the  3ame  approximations  as  before  are  made,  the  average 
absorption  cross  section  for  energy  group  i  becomes 


1,.  -  tV/^i 

^j£s}1  avg 


Vi/f2sJi  avg 


III*   Spatial  Integration  of  the  Group  Diffusion  Equations. 
Equation  (1.5)  describes  the  neutron  behavior  in 
energy  group  1*   Substitution  in  (1*5)  of  (2,5)  plus  the 
following  relationships  developed  previously 

qi0UT  -    Uli^  {u6) 


Dl  = 


(  1        ) 

(TpTT?i  avg 


NV1 

(q)lavg:T7TXI 


(20) 


(2.'+a) 


U 


l(/^i   avg 


gives 
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-n^N^+jr.i+5-qji, 


(J^sJl  avg' 


)  NV1 


-   OUT   , 


ys(r,u 


)du 


(3-D 


For  spherical  geometry,  the  following  substitutions  will 

2 

simplify  (3»l)  by  eliminating  the  Laplacian  operator  V  : 

1*'n  =  rnNVl  (3. 2a) 


1  OUT  -    OUT 
Si  ,.  -  rnqi 

^nr  rn/3<rn*u)du 


(3.2b) 
(3.2c) 


A  new  notation  is  adopted.  Let  i  refer  to  the  i^th 
energy  group  as  previously  and  the  subscript  n  refer  to  the 
n'^h  space  point  located  at  a  distance  r  from  the  origin  of 
the  spherical  coordinates.   As  before,  the  source  term  ^ 
represents  neutrons  entering  group  i  at  radius  rn  as  a 
result  of  inelastic  scattering  and  fissions  in  all  groups. 
Use  (3.2)  in  (3.1)  to  get 


-% 


.2  ip 


n 


+ 


dr 


^-a^ 


+ 


L 


f^  (_1) 

^/«£s'l  avg 


^n  =^  0UT 


^  *n 


(3.3) 


The  space  points  are  chosen  to  divide  the  spatial  region 
into  a  3et  of  subintervals  of  width  ^r.   The  finite  dif- 
ference approximation  . 


<*2  fa     =^L±^Vi 

K   dr2  ^rn         (Ar)2 


-  2  1F 


H 


(3.4) 


and  Equation  (3*3)  gives  the  recurrence  relation 


where 
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lrn  +  l  =  *l\-tF».l-\  (3.5) 


bl  =   la,  +  f  T7    1    x    )  <3.5b) 


1-l^OUT     +  19J  (3.;o) 


Equation  (3.5)  can  be  verified  by  direct  substitution  of 
Equation  (3.^)  into  Equation  (3.3).   The  above  substitutions 
follow  a  method  developed  by  Ehrlich  and  Hurwitz. 

a.  Boundary  Conditions*  The  boundary  conditions 

for  the  numerical  solution  of  Equation  (3.5)  are  that 
*F  s  0  at  r  —  0  and  also  at  the  extrapolated  end  point  out- 
side the  assembly.  Another  condition  is  that  the  neutron 
flux  and  current  density  must  be  made  continuous  at 
boundaries  between  regions  of  different  material  properties* 
Suitable  choices  of  /\  r  must  be  made  for  each  region  so  that 
the  boundaries  will  occur  at  space  points.  At  a  boundary 
point,  M,  located  at  radius  r  -  rjj  between  two  different 
media $   these  boundary  conditions  give  (primes  refer  to  the 
outer  region)  .. 

iFj4  -  1fV'  for  continuity  of  flux, 

^M  =  ^*M»  for  continuity  of  current  density, 
where 
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j  =  -D  &m  z  . jjJ-P  +  r  &> 


dr 


dr' 


It  oan  be   shown  that      (see  Appendix  C) 


(-F  +  r  g) 


^  - 


V£"  tFMMr  -  FK  -   Ar)   -  I  FM       <3.5d) 


2^r 

~  6  ^FM-h^r~'"  FM  -  A  r^ 
Therefore,  the  boundary  condition  for  neutren  current 

density,  J  s  J1 ,  Is 
f  FM  +  £  (FM^r  +  F„  .^r)  -Jk-  (FM^r  -  FM  ^r) 


-  D» 


3 

:f« 


f  F'M  +  £<FWr'+FWr«>   -«3U 


2Ar' 


-  F      ) 

M  +Ar    M-  Ar 


(3.6) 


b.  Method  of  Solution*  The  recurrence  relation 
for  ^Fn,  Equation  (3«5)»  will  be  solved  using  a  procedure 
suggested  by  Ehrllch  and  Hurwitz.   This  procedure  is 
particularly  convenient  for  use  v/ith  the  WWI  digital  computer. 

To  solve  Equation  (3.5)  subject  to  the  above 
boundary  conditions,  new  variables  c*n   and  0  n  will  be  intro- 
duced. These  variables  are  related  to  the  homogeneous 
(in  s  0)  solution  of  (3.5) »  An,  by  the  defining  relation- 
ships 

■  &±i-  and  9      Z      f"     A«  ll 


\ 


c^ 


When  the  above  expression  for  oC   is  substituted  into  (3*5) > 
with  In  =  0,  It  is  found  that 

(3.7) 


n  + 1  :  *  sen 
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Th  e  recurrence   relationship  for  the/3's   i 


B 


which  can  be  verified  by  direct  substitution  of  the  definition 
of  0 n   into  (3.8).  A  combination  of  Equations  (3.5),  (3.7), 
and  (3«8)  gives  the  following  relationship 

Fn  =  Jn^  +  ""  (3-9) 

^n+l 

This  is  proved  as  follows: 

Substitute  (3.7)  and  (3.8)  for  <X  and  &   into  (3.9) 
and  combine  terras.  This  gives 


^  |Fn  +  #„-! 


=Fnk  "  Fn+1  "   I„ 


3y  Equation  (3.5),  the  right  side  becomes  F  -,. 
This  is  then  Equation  (3.9),  when  n  is  replaced  by  n+1. 

Equation  (3.7)  is  used  to  solve  for  the  <X's, 
rtlng  with  °^2  =  &t  since  A0  =  0  and  therefore, from  the 
definition  of  oe  f     c/  -CD.    The  next  step  is  to  use  (3*3)  to 
solve  for  the  /9  's,  starting  with  (B     =  I-.  Finally  (3.9) 
is  solved  for  the  F  's  from  the  outside  to  the  center  of  the 
spherical  reactor  assembly,  starting  with  F»  :  C  at  the 
extrapolated  end  point. 

The  ecuations  for  ex  and  f3    have  the  same  mathe- 
matical form  but  different  constants  in  regions  of  different 
material  properties.  At  the  boundary  between  regions,  o^  and 
0    are  single  valued  in  order  to  maintain  continuity  of  flux. 
Thus  Equation  (3.9)  is  unchanged  at  a  boundary.  At  the 
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boundary  point  M,  at  radius  r.  ,  separating  two  different 
media,  the  boundary  condition  for  continuity  of  current 
density  can  be  satisfied  by  using  Equation  (3.6).  The  sub- 
stitution of  the  recurrence  relationships  for  o<  and  /5  into 
(3.6)  leads  to  the  following  relationships  across  boundary 
point  K  (primes  refer  to  the  outer  media): 


£r*.    1 


(3.10) 


(3.11) 


where 


s  = 

D 

P  ■  -*-(-¥^  -  £H  =  -f-  ^f       '  (3.12) 


A  =  pk  +   1  + 


d   .  l\2Ar' 


c.  Thermal  G-roup.  The  thermal  neutrons  are  as- 
sured to  be  in  a  group  of  essentially  constant  lethargy.  It 
is  assumed  that  neutrons  enter  the  thermal  group  only  by 
elastic  scattering  from  the  group  of  next  higher  energy. 
The  thermal  group  is  treated  the  sa:::e  as  other  groups  except 
for  changing  the  definition  of  KVm  (where  the  subscript  T 
stands  for  the  thermal  group)  •  l(1or  the  thermal  group 


*  for  the  case  considered  here  with  15  space  points  in 
the  core  so  that  A  r  =  R/15,  where  R  is  the  core  radius 
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NVT  s  fa 

&•   Energy  Gr o\xr>a *  Energy  Group  I  is  chosen  from 
10  Mev  to  0,1  Mev  (u  -  0  to  u  s  4.606).   Glasstone^  plots 
the  fission  neutron  energy  spectrum ,  and  it  can  be  seen  from 
this  plot  that  very  little  error  is  introduced  by  assuming 
that  all  neutrons  released  in  fission  are  born  into  G-roup  I* 
Groups  II  and  III  are  chosen  to  give  approximately  equal 
lethargy  widths  down  to  a  thermal  energy  corresponding  to 
1000°F. 

The  energy  groups  will  be  as  follows? 
2£2ffi    jLethargy,  u_       Energy,  E_     Lethargy  Width 
I       0  *  b.6o6  10^0,1  Mev         4.606 

II    4.606  -  10.81       100-0.2  Kev  6.204 

III   10.81  -  18.80       200-0,070  ev       .  7.990 
T       18.80  0.070  ev  (1000°F)       o 

e.  Inelastic  Spattering*  Provision  is  made  for 
considering  inelastic  scattering  into  an  energy  group  by  the 


Js(r,u) 


definition  of  the  term  J   S(r*u)du  in  the  age  diffusion 
equation.   This  term  is  defined  to  be  composed  of  fission 
neutrons  born  into  group  i  plus  those  inelastically 
•cattered  into  group  i  from  higher  energies.   Since  a 
neutron  must  possess  a  minimum  of  about  0.1  Mev  of  energy  if 
it  is  to  be  Involved  in  an  inelastic  scattering  process, 
Inelastic  scattering  occurs  only  in  Group  I.   The  group 
widths  are  large  so  there  will  be  no  inelastic  scattering 
across  groups.  This  can  be  seen  from  plots  of  the  energy 
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dlstribution  of  inelastically  scattered  neutrons  as  given  in 

the  Reactor  Handbook,  Reference  (?)• 

IV,  Programming*  for  the  Whirl vl nd  £  Cpnputer* 

Whirlwind  I  (WWI)  Is  a  digital  computer  capable 
of  performing  simole  mathematical  operations  (addition,  sub- 
traction, multiplication,  and  dlrlslon)  with  extreme 
rapidity  under  the  control  of  a  set  of  instructions,  called 
the  program,  which  is  stored  in  the  computer.   A  detailed 
account  of  WWI  computer  operation  is  given  in  Chapter  II, 
Section  II,  of  this  thesis.   For  additional  information, 
the  reader  is  referred  to  Reference  (l) ,  which  completely 
covers  WWI  operating  techniques, 

a»  Numerical  Operations*  3ection  III  of  this 
chapter  develops  the  relationships  which  are  used  in 
programming  the  WWI  solution.   The  steps  of  the  solution  are 
outlined  below: 

(1)  Start  by  assuming  the  spatial  distribution 

of  fission  neutrons  born  into  Group  I  as  a  result  of  fissions 
in  all  groups.   (This  assumed  distribution  will  be  Improved 
to  the  correct  value  by  successive  Iterations.) 

(2)  From  the  assumed  spatial  distribution  of  fis- 
sion neutrons  into  Group  I,  which  la  J   3(rn#u)du,  use  (3.2c) 
to  get  3n  i 

Isn  -  rn  J  S(rn»u)du. 

(3)  Calculate   the   following  constants   for  Group  If 
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using  Equations  (3.5)  and  (3.12): 

bi  =  Zai  + 


'>"'tH 


S  "   D 


avg 

D' 


P  ■ 


-   1  /Ar1   Ar' 


2s  (Ar    3r 


rJ 


MAr) 

k  =  2  -f-  •  - 


2 


Di 


T  I 

(4)  Calculate  ©<n  starting  with  oC^  ■  k  and  using 

Equation  (3-7)  ,  c<  -,  ■  k  -  — t~-   The  value  of  k  must  be 
'   n*l       oCn 

correct  for  the  region  considered.  At  the  boundary  between 
core  and  reflector,  where  n  =  15,  boundary  condition  (3*10) 
is  used,  c*'   =  K TC-'—T^-- 

(5)  Calculate  I/0n  starting  with  fB  ^  =   1^  and 
usi       tion  (3»8)i  ^T,  =  4"  *r>j  where 

III   =  0  for  16  -  n  ^  25 

Ij. 

At  the  boundary,  n  =  15,  Equation  (3.11)  applies, 

Ar*    #14 


*« "  wis  +  (i  *  ^  + 


sAr  exf,,- 


T  I 

(6)  Solve  for  x?      starting  with  F25  =  °  and 
working  toward  the  center  of  the  spherical  reactor  assembly 
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using  Equation  (3»9)> 

(7)  Calculate  the  Inelastic  scattering,   3n  -  , 
from  Group  I  into  Group  II  by  use  of  the  relationship 

in.  x 

where  3C  is  the  probability  that  a  neutron  scattered  in- 

elastically  in  Group  I  will  enter  Group  II • 

(8)  Repeat  successively  steps  (2)  through  (6)  for 
Groups  II,  III,  and  T.   Note  that  there  are  two  changes  to 

the  procedure  used  with  Group  I:   First,  the  input  term, 

II 

3  ,  is  composed  only  of  lnelastically  scattered  neutrons 

entering  Group  II  (since  all  fission  neutrons  are  born  in 

Group  I),  and  $n   is  zero  for  Groups  III  and  T;  second,  the 

calculation  of  xIn  by  Equation  (3«5o)> 

2 


iin  .  L&xl 

*n  - 


Dl 


*n     n 


i— 1  OUT 
includes  elastic  scattering,  n     ,  out  of  the  group  of 

next  higher  energy. 

(9)   When  Fn  ha3  "been  calculated  for  all  energy 

I 
groups,  calculate  a  new  source  term,  5n,  of  fission  neutrons 

for  Group  I  from  the  relationship 

-SnZ  2Z    UZf>!  avg  ipn' 
1=1 

where  (yZjO^  av?.  is  the  average  value  in  group  !•   Thus, 
*Sn  Is  proportional  to  the  fission  power  distribution  in  the 
reactoro 
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(10)  The  new  input  terra,   Sn,  is  normalized  and 
used  as  the  basis  for  a  new  iteration,  which  repeats  the 
above  procedure  starting  with  step  (3). 

b.  Iteration  Procedure.  The  iteration  procedure 
has  been  divided  into  two  steps:  (1)  relaxing  from  the 
assumed  input  of  fission  neutrons,   Sn,  to  the  correct 
spatial  distribution;  and  (2)  changing  the  radius,  R,  of 
the  new  reactor  core  to  obtain  criticallty. 

Let  the  superscript  m  indicate  the  number  of  the 
iteration,  with  m  =  0  indicating  initially  assumed  values. 
Let  a  circumflex  indicate  values  which  have  been  normalized. 

Step  (1)  corrects  the  spatial  distribution  of  Sn 
as  follows: 

T    m 

All  iterated  input  terms  Sn  into  Group  I,  starting 

from  the  assumed  Sn°,  are  normalized  with  respect  to  an 

*  I   m-l 
arbitrarily  chosen  core  spatial  point,  N,  where*  SN    is 

normalized  to  unity.  A  normalizing  factor  %   is  defined  such 

that 

l£  m-l 

In  m      1q  m 


Lo  III         J-o 

5N        faW 


and 


n       u   n 


where  N  =  the  chosen  spatial  point  for  normalization  to  unity 
m  =  the  number  of  the  iteration 
Snra  =  the  calculated  input  term  to  G-roup  I  at  spatial 
point  n 
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1A  m  - 
3  '  -  normalized  source  team 
n 

ja  m— 1  — 
3„    -  l  s  source  term  of  previous  Iteration  normalized 

to  unity  at  space  point  N 

To  determine  when  to  stop  iterating,  the  normalized 

Ia  m 
Input  term,  3   ,  is  compared  at  each  spatial  point,  n,  with 


the  preceding  iteration  (see  figure  below).   It  is 


Ism 

n 


normalization  ^ 

point      I 


3 


.  13  15 


a  m— 1  a  m 

specified  that  at  eaoh  spatial  point  the  ratio  3n   /3n 

mu3t  differ  from  unity  by  no  more  than  a  small  value,  £  , 
to  permit  the  program  to  proceed  to  the  next  step.   Whenever 
the  value  £  Is  exceeded  at  any  space  point,  the  program  T-rill 
reiterate.   This  is  equivalent  to  requiring  that  the  solution 
for  the  core  power  distribution  must  have  converged  suffi- 
ciently to  permit  the  solution  from  one  iteration  to  lie  en- 
tirely within  an  envelope  determined  by  £    around  the  so- 
lution of  the  preceding  iteration. 

The  criterion  for  choosing  between  iterating  or 
continuing  to  the  next  step  is  based  upon  the  ability  of  >/WI 
computer  to  choose  one  course  of  action  for  a  positive 
number  and  a  different  one  for  a  negative  number  (conditional 
control  transfer).   Provision  is  made  for  the  criterion  to 
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ID-1  / 


operate,  Irrespective  of  whether  3n   /3n  is  greater  or 
less  than  unity,  by  the  following  schemes 

(a)  When  S^ is  greater  than  unity: 

8  m 
n 

1  -  n    -/-£  s  0,  not  acceptable  (reiterate) 

sn 

/\ 

S  nj-1 

1  -  n    )  £  z*0*  acceptable  (prooeed  to  the 
/\  m 

3n  next  step) 

3  nwl 

(b)  When  -2 is  less  than  unity: 

s  m 

n 

s  m"1 

-  n    ~1  -f-  €  <C  0,  not  acceptable  (reiterate) 

£  m 

sn 

°n    -1  +6  /  0,  acceptable  (proceed; 
*   m 

Step  (2)  of  the  iteration  procedure  involves  changing  the 
core  radius  to  get  criticallty<>   The  procedure  is  as  follows: 

Step  (l)  above  has  required  that  a  plot  of  the 

^  m 
innut  term  Sn  must  reach  a  stable  shape  before  procedlng  to 

step  (2).   In  step  (2)  the  last  calculated  value  of  ^  k>Y 

(Equation  ^.1)  is  used  to  test  for  over  or  under  criticality. 

A  ft  less  than  unity  means  that  the  number  of  fission 

neutrons  entering  energy  G-roup  I  increases  from  iteration  to 

iteration  and,  therefore,  that  the  multiplication  factor,  k, 
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ia  greater  than  unity  (and  vice  versa).   For  the  critical 
condition,  k  -  1,  the^source  term  remains  constant  for  suc- 
cessive iterations,  and  #   therefore  is  unity. 

If  <y  is  not  within  a  small  value  €    of  unity,  the 
program  goes  through  another  iteration  using  a  new  value  of 
the  core  radius,  R,  chosen  to  come  closer  to  criticality. 
The  new  radius  for  iteration  m  is 

Rm  Z    ^m-1  Rm-1 
Thus,  the  core  radius  is  increased  if  k  <C  1 ,  and  vice  versa. 
The  criteria  for  determining  whether  the  radius  is  correct 
or  not  is: 

(a)  When  ^  is  greater  than  unity: 

6   -  #+  1^0     not  acceptable  (modify  radius 

and  reiterate) 

6  -^+  1  /*  0     acceptable  (proceed  to  the 

next  step) 

(b)  When  Y  is  less  than  unity: 

y-  1 -/■  6  <0     not  acceptable 
y~   1  +■  6  }  0     acceptable 
When  the  iteration  criteria  of  both  step  (l)  and 
step  (2)  are  satisfied,  the  program  proceeds  to  the  final 
step  of  calculating  and  printing  the  core  power  distribution, 
the  critical  radius,  and  the  neutron  flux  distribution  in 
each  energy  group. 
v*  Hoy  Sheet . 

The  program  for  the  solution  of  the  reactor  problem 
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has  been  described  in  3ection  IV.  A  flow  sheet  will  be 
developed  to  show  schematically  the  operations  discussed 
in  3ection  IV.   The  following  conventions  will  be  used: 

A  square  box  standing  on  one 
edge  denotes  a  subroutine  for 
calculating  the  quantity  within 
the  box. 


Circles  enclosing  numbers  are 
used  to  specify  corresponding 
parts  of  the  flow  sheets.   The 
numbers  are  called  "addresses. M 
The  circles  are  called  connectors 
and  are  put  at  those  places  In  a 
flow  diagram  to  which  reference 
will  be  made  at  some  other  point. 


A  rectangular  box  denotes  that  the 
quantity  Inside  is  to  be  calcu- 
lated or  the  operations  within 
are  to  be  performed.   Frequently 
the  formula  for  calculating  the 
quantity  is  given. 


^V> 


A  :  b 


j 


a 


This  symbol  Indicates  a  conditional 
transfer  to  a  If  A  >  B  and  to  b 
If  A  —  B. 
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Subscrlpts  and  superscripts  will  have  the  follow- 
ing definitions: 

i  denotes  an  energy  group* 

n  denotes  a  space  point  at  radius  rn« 

m  denotes  the  number  of  an  iteration,  where  m  r  0* 

denotes  values  initially  assumed* 
N  denotes  an  arbitrary  space  point,  n  -  N,  chosen 

as  a  basis  for  normalization  procedures. 

■>  the  symbol  a«Vb  means  that  the  quantity  a  replaces 
the  quantity  b  in  the  storage  register  in  which 
b  is  stored. 
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FLOW  SHEET  FOR  FOUR  GROUPS 
Main  Program 


V- 


Stored 
Constants 


i  =  l 


Put  Gr.  i 
constants  into 
operating  registers 


<D 


<EK^y<^}^<^) 


> 


II 


Snin  =  Ui(xIin.)IIFn 


n+1 


JL-  starting  at  F '  =  0  at  reflector  outer 

boundary . 


FLOW    SHE!  T   FOl  *ROoFS 

\J1  -Fro  jrar.    (Continued) 


^— 


n 


-  +6  :0 


> 


:-l 


n 


s- 

"*n 

> 

> 

t. 

' 

*■» 

a  - 

n 

— >■  o 

n 

1=1 

1  +f  :0 

-. / 


to 


© £    -y+g'-y^ (^   1   -hf:o)^-(7)— 


Print  <T ,    R 
R"'+1    =  JmRn  -> 


A 

o 


A     -1 

^    n 


1=1 


0-© 


!  Print 
I  (all   n) 


s-l 


(S^C^nV©- 


-^ 


irmt 

Ml 

I       in) 


'm 


n 
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FLOVJ    SiiEEl 

Grouo  Constants  Subrout: 


©■ 


3   sR 


k  =  2  4  _M^' 


ts> 


£> 


p*^*^K  +(i--i-j%- 


© 


(8> 


Q(n  Subroutine 


« =  ,. 


..   .  ,J  1 ■  .    ■■       ■■  i    .ii       ■  -    ■  ,-i.-,  —.,■—'-— j  y 


n  4    1    -* 


©■ 


C<  =     X     - 

16 


Ail    i 


iAr      °<I 


1=17 -(Tv 


1 


<*n-l 


-<El>Hi) 


boundary 
condition 


+  1    -*•  n 
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FLO'rt    SHEET   J?"Oirt  FOUR   OKOUPS 


(3    Subroutine 
n 


Q 


O 


in  =  g2i,  fen 


13 


o< 


■■  +  I 


n 


n 


<n7TZ>_>-<b) 


n    +   1    ->•  n 


o 


&- 


l^Ar'  .  ^lA. 


15  -    ^)2>_   •  ft5  =  >h!  +    ^  S 


.i=   16 


Kf> 


boundary  condition 
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FLOtt    GIIEET  FOR  FOUR   GROUPS 


Pn  Bui 


broutine 


e-E£ 


7^->_® 


0H 


I15.  I{5,^5 .  PI15  +fi-  +$E>i 


+ 


l'4r  ' 


.Ji 


1  A 


£H 


boandary  condition 


a 


On-l 


n 


c* 


+       I 


> 


nfl 


© 


ex 


KD 


O 


.is  forr.  of  I  comes  from  combination  of  Equations  (1.6), 
(2.4a),  (3.2),  and  (3.5c)  wither  *  R/l5,  wh    i  =  core 
radius.  (There  are  15  core  space  points.) 
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flow  s::fft  for  four  groups 


L0  Subroutine 


(oY4n=TT-(T) — 0n  =  ZDl 


^-^^ ^^  n   R  +   (n-1 


(n-15)CAr*) 


<5^-<T) 


fn+1  *n 


© 


Flad 

Pi 
p2 
bl 


11 
cl 

c2 
c3 


Memory  Assignment  (See  W'.VI  program  following) 
Content 9 

Storage  of  numbers  used  in' numerical  operations 
Temporary  storage 
R,  (v2f)IIf  (j/£f)ni,  (^^f)T*  T»  <Ar') 

(X*ln>i.  CX^1,  6  ,  (y^f)I 

Used  to  indicate  number  of  iterations 

Group  I  constants:   D,  D« ,  f,  f»,  (j^s)L   avg» 

'f^iavg.  ^a'  V'  ^i 

Group  II  constants:   (Same) 
Group  III  constants:   (Same) 


*  A  FLAD  is  a  floating  address  used  to  locate  a  particular 
set  of  storage  registers. 
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Flad  Contents 

cb  Group  T  constants:   (Same) 

I  o 
si  Assumed  source   S 

s2  Calculated  source  term  from  fissions  in  all 

T 

groups   2  (*^fFn)     "   3n 
3-z1 

fl  JFn       l-n^25 

f2  nFn     l^n^25 

III 

O  Fni£n^25 

fk  TF  l£n*25 

n 

zl  Operating  register  for  storage  of  registers  cl 

when  group  i  is  being  used 

z2  Temporary  storage  for  o< 

n 

1 

z3  Temporary  storage  for  0 

n 

zz*  Temporary  storage  for    Fn 

xl  Temporary  storage  for  inelastic  scatt,   ^bn  in 

rl  Group  constants  subroutine 

r2  «=*  subroutine 

n 

r3  1/3    subroutine 

n 

r^  &n   subroutine 

r5  ^-Fn  calculation  stored  in  fl 

r6  Fn  calculation  stored  in  f2 

ITT 
r7  Fn  calculation  stored  in  f} 

T 

r8  Fn  calculation  stored  in  f^ 

r9,  rlO        Calculates  (  3n^ln.  inelastic  scattering  and 

store.9  in  xl 
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Flad  Contents 

rl2  10n  subroutine 

frl3,  810       Iterates  for  correct  spatial  distribution  of 

flux 

yl4,  all       Changes  radius  and  reiterates  until  criticality 

is  achieved  (both  rl3  and  rl^  will  perform  a 
maximum  of  50  iterations  before  printing  out 
the  answer  so  as  to  avoid  a  waste  of  computer 

time  should  the  solution  be  divergent) 

I  o 
x2  Calculates   Sn  from  the  assumed  Sn  where 

I3°  =  rn(3n) 

Control  for  Group  I  calculations 

Control  for  Group  II  calculations 

Control  for  Group  III  calculations 

Control  for  Group  T  calculations 


al  and 

yi 

a2  and 

y2 

a 3  and 

73 

a 4  and 

y4 

s8 

39 

slO 

811 

s!2 

Calculates  ^-  1  where  N  is  spree  point 


1C 


q  m 
3N 


Normalize  Sn  -  o  3n 
(See  rl3) 
(See  rl4) 

Puts  0n  in  r5>  r6,  r7,  r8  in  place  of  Fn» 
Prints  core  power  distribution,  flux  distri- 
butions, critical  radius,  o   ,  and  number  of 
Iterations 
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WHIRLVJIND  PROGRAM 


fc  TAPE  374-223-008  A  HOOVER 
(24,6) 


t=p2 

pi, +1.0 

+7.0 

4-13.0 

+50.0 

+0.0 

+20. 

11, +0.0 

+1.039 

+0.00rl 

+0.522 

+0.0 

+1.05 

+0.12484 
+0.03784 

si, +1.33 
♦1.15 

+0.93 

s2,+0.0 

fl,+0.0 

f3,+o.o 

Si, +0.0 

z3»+o.o 

xi,+o.o 

idvzl+4 

itszl+20 

idvbl 

lmrzl+18 

lmrzi+10 

ladpl+2 

lmrzl+28 

icspl 

idvbl 

icazl+24 

iadzl+24 

Idvbl 

IscO 

itsz2+30+c 

icabl 

imrsl 


+2.0 

+8.0 

♦14.0 

p2,+0.0 

bl,+60.0 

+2.0 

+0.0 

+0.G21 

+1.528 
+0.16119 

+0.03813 

c4, +1.190 
+0.1282 

+1.31 
+1.10 
+1.14 

DITTO 

DITTO 

DITTO 

DITTO 

DITTO 

DITTO 

irarzl+S 

icapl+20 

itszi+22 

idvzl 

iadzl+l4 

it sz 1+28 

idvpl+2 

idvzl+20 

iad2t 

itsz2 

Itsz2+2+c 

idvzl+20 

icr9 

ldvpl+28 


||  itsz3 
Its2V     icaz3+c 
lctllr    icazl+22 
lmrbl+10  idvbl 
ltsz3+2o  isc5 
itsz3+30+c 
icabl     idvpi+28 
icazl+36  imrzl+4o 
imrt      itsz3 


48r 
48r 
40r 
48r 
48r 


+3.0 
+9.0 
+15.0 
+0.0 

+0.001073 
+0.0 

+0.0 

+0.68 
+0.0 
+1.08 
,♦1.184 

+0.12968 

+0.0354 

+1.29 

+1.05 

+i.45 

28rj  +0.0 
+0.0 
+0.0 
+0.0 
+0.0 
+0.0 

ladzi+12 

imrpl+2 

icabl 

iadpl+2 

itGZi+26 

icabl+iO 

ltat 

iadpl 

iadt 

isc4 

idvz2+26 
icspl 
ict20r 
Imrbl 

isc4 

idvz2+c 

irart 

idvzl+20 

icrlO 

lct33r 
irarbl 
idvzl+32 
isc4 


+4.0 

+10.0 

+16.0 

+0.0 

+0.0 144^ 

+0.00 

+0.0 
+0.02158 
+0.2171 
+0.06535 

+0.008951 
+0.544 

+0.003914 

+1.26 

+1.00 


+5.0 

+11.0 

+0.0 

+0.0 
+0.052472 

+0.03 

cl,+2.209 
+0.06884 

c  ,+0.5756 
+0.13231  +0.000523 


+6.0 

+12.0 

+5C0 

+0.0 
+0.000219 


+C.537 
+0.000522 

+1.0 

+1.22 
+0.95 


f2,+0.0 

f4,+0.0 

z2,+0.0 

z4,+0.0 

rl,ita54r 

itszl+i8 

idvpl+4 

idvpl+28 

itszl+24 

imrbl+10 

idvpl+4 

Icazl+22 

lmrpl+2 

itszl+30 

icrl3 

ict5r 

imrbl+10 

idvz2+28+c 

ispO 

idvpl+28 


DITTO 
DITTO 
DITTO 
DITTO 

icazl+2 

imrbl+10 

imrbl 

icazl+l6 

imrbl+10 

idvbl 

imrzl+24 

idvpi+4 

ispO 

icspl 

icspl 

Iadzl+30 

r3,ita37r 
idvzl 


icri3  icat 

iad2t  itsz3+2+c 

imrsl+28  its2t 

imrz3+26  idvz2+26 
icaz3+28+c 

ispO  r4,ita75r 

idvpl+28  idvzl 

its2t  imrz4 

icr!3  ica2t 


+1.66 

+1.0 
+0.0 
+1.20 
+O.92 


48r|  +0.0 

48r  +0.0 

46r)  +0.0 

48r|  +0.0 

icazl+16 

idvzl 

idvzl+20 

idvpl+28 

idvzl+6 

idvzl+2 

iudpl  • 

its2t 

imrbl+10 

r2,ita25r 

Idvz2+c 

imrpl+28 

itsz2+28 

iadzl+28 

itst 
imrsl+2+c 

leap 1+28 

iad2t 

idvz2+28+c 

it3t 

iadxl 

lmrz4+2+c 
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VJhirlwind 

Program 

iadxl+2+c 

irurt 

its4t 

icaz3+c 

idvz2+c 

iad4t 

itsz3+2+c 

ictl7r 

ica2t 

imrz4H-28 

iadxl+28 

imrt 

lmrzl+22 

its4t 

icabl+10 

imrbl+10 

idvzl+2 

itst 

icazl+38 

iinrzl+4o 

idvzi+34 

its2t 

imrz4+28 

iadxl+28 

imrt 

lts6t 

icabl+10 

ldvpl+10 

idvbl 

itsBt 

leap! 

idvpl+2 

iad3t 

imr6t 

its6t 

icabl+10 

lmrpl+28 

idvbl 

idvzl+20 

imrz3+26 

idvz2+26 

iad4t 

iadot 

itsz3+28 

i3c5 

icrlO 

ica2t 

imrz4+30+< 

c 

iadxl+30+c 

imrt 

its8t 

icaz3+28+c 

idvz2+2C+c 

iad8t 

itsz3+30+c 

ict65r 

ispO 

r5,itallr 

icapl+32 

itsfl+48 

iscG 

icr24 

icafl+43+ 

c 

iadz3+46+ 

c 

idvz2+46+c 

itsfl+46+. 

o 

icd2 

ict5r 

ispO 

lllillill 

icr24 

r6,itallr 

lcapl+32 

itsf2+48 

isc6 

leaf  3*48+0 

iads3+46+c 

idvs2+46fc 

itsf2+4G+i 

T 

icd2 

ict5r 

ispO 

r7,itallr 

leap 1+3 2 

"'itofS+^S 

isc6 

icr24 

icaf3-i 

c 

iadz3+46+c 

idvz2+46+ 

c 

ltsf3+46+ 

c 

icd2 

lot5r 

ispO 

r8,itallr 

icapl+32 

itsf4+48 

isc6 

icr24 

icaf4+48+ 

c 

iadz3+46+c 

idvz2+46+< 

1 

itsf4+46+ 

c 

icd2 

ict3r 

IspO 

r9,ita!3r 

isc3 

icrl5 

icafl+c 

imrbl+12 

itsxi+c 

ict3r 

isc5 

icrlO 

icafl+30+c 

imrb  14-14 

itsxl+30+ 

c 

ict9r 

ispO 

rlO,ita6r 

lsc2 

lcr25 

icapl+32 

itsxl+c 

ict3r 

ispO 

rl2,ita2i; 

p 

icapl+28 

idvbl 

itst 

isc3      icr!5 

icaz4+c 

imrt 

idvpl+c 

itsz4+c 

ict6r 

isc5 

icrlO 

icabl+10 

irnrpl+c 

iadbl 

its2t 

icaz4+30+c 

ldv2t 

itsz4+30+ 

c 

ictl3r 

ispO 

rl3#isc3 

lcrl5 

icas2+c 

itssl+c 

ict2r 

icail 

iadpl 

itsil 

icail+4 
rl4,isc3 

lad  pi 

itsil+4 

i  sup  1+3'^ 

icpal 

ispsil 

lcrl5 

icas2+c 

itssl+c 

ict2r 

icat 

iMOA-123.1234t 

isupl 

idvpl 

iadpl 

imrbl 

itsbl 

1MOA- 123. 1234c 

leal 1+2 

iadpl 

itsil+2 

icaii+6 

iadpl 

itsii+6 

i sup 1+3 6 

icpal 

ispsl2 

x2#isc3 

icrl5 

icasl+c 

ldvpl+18 

imrpi+c 

itssl+c 

ict2r 

lillflllll 

111 

al,iscO 

icr9 

icacl+c 

itszl+c 

ictal+2 

isprl 

lspr2 

ispr3 

ispr5 

ispr9 

yl,isc3 

icrl5 

1 cab 1+18 

imrf 1+c 

itss2+c 

ictyl+2 

a2#iscO 

icr'9 

icac2+c 

itszl+c 

icta2+2 

iscl 

lcr4 

icacl+4+c 

ltszl+32+c 

i eta 2+7 

icacl+16 

itszl+4o 

lsc2 

icr25 

icafi+c 

itsz4+c 

icta2+l4 

ispri 

ispr2 

ispr4 

lspr6 

isprlO 

y2,isc3 

icrl5 

i cab 1+2 

imrf2+c 

iads2+c 

itss2+c 

icty2+2 

a3,iscO 

icr9 

lcac3+c 

itszl+c 

icta3+2 

iscl 

icr4 

icac2+4+c 

ltszl+32+c 

icta3+7 
itsz4+c 

icac2+l6 

itszl+4o 

isc2 

lcr25 
ispr4 

icaf2+c 

icta3+l4 

isprl 

ispr2 

lspr7 

y3#isc3 

icr!5 

icabl+4 

irnrf3+c 
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hirlwlnd 

Pre 

iads2+c 

3132+C 

by3+2 

,iscO 

ler9 

lcac4+c 

ituzl+c 

icta4+2 

iSCl 

icr4 

icac3-fi+«»-c 

itszl+32+c 

icta4+7 

3cac^+l6 

icta^i; 

itszl+^O 

ice  2 

Icr25 

icaf3+c 

itsz^+c 

Isprl 

lspr2 

ispr4 

isprB 

y4,l3C3 

Icrl5 

1 cab 1+6 

ImrfJf+c 

:G2H-C 

ltB»2+C 

lcty4+2 

s  jicapi 

idvs2+- 

itst 

,Isc3 

icrijj 

icat 

irar3  2+c 

it3G2+c 

Icts9+2 

siC,isc3 

icr!5 

icasi+c 

idvs2+c 

its2t 

icapl 

iadbl+16 

Isu2t 

Icprl3 

icspl 

iadbl+16 

iad2t 

icpr. 
Icpri4 

ictslO+2 

■j.#icabl+ 

I  out 

iadpi 

tat 

pi 

iadbl+' 

icprl1* 

IFOR 

+c 

+25 

,:.: 

lcr2p 

icafl+c 

itaz^  ; 

ict2r 

lcprl2 

is 

lcaz4+c 

Jfl+C 

.  Ti' 

icc2 

Icaf2+c 

itszW 

ictllr 

isprl2 

1dc2 

icr25 

icaz^+c 

ifcsf2+e 

ictl/r 

1gc2 

icaf3+c 

itsz^+e 

ict21 

laprl2 

Isc2 

Icaz^+c 

itsf3+c 

lct2( 

isc2 

icaf^+c 

ltGZ2*4-C 

ict30r 

ispr 

isc2 

icaz4+c 

ItafHc 

ict35r- 

iBC 

icafi+o 

idv 

itsfi-fc 

1MOA-P123 

.1234; 

ict^Or 

isc2 

icaf2+c 

idv 

+c 

iMOA 

lct; 

1gc2 

IcafJ+c 

Idv 

+c 

ict£2r 
|4t 

icaf4 

Id  ! 

leaf  ^-t-c 

IdvH 

Itsf4+C 

IMOA 

icaii 

DA 

icaii+2 

IMOA 

icabl 

DA 

isc3 

icr!5 

icarj2+c 

BPpI+l8 

idvpi+o 

b73r 

- 

OP 

.  :T  AJ 

.*.  j-i— 
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Table  IV  of  Appendix  B  gives  the  Whirlwind  input  data 
for  runs  (1)  through  (10).  The  Whirlwind  program  for  run  (8) 
has  been  given  on  the  proceeding  pages.  The  storage  assignments 
for  the  input  data  are  as  follows  (using  run  (8)  as  an 
example) : 


Numerical 
i'LAD       Value 


bl, 


cl, 


c2, 

o3, 

cA, 


60.0 

0.001073 
0.01444 
0.052472 
20.0 
2.0 
6.0 
0.0 

0.03 
0.000219 

2.209 
1.039 

0.621 
0.68 
0.02158 
0.06684 

0. 000106 

0.0 

0.2171 


Description  of  quantity 


R,  assumed  core  radius  (cm) 

Y\t  in 

U£f  )T 

T,  reflector  thickness  (cm) 
^r*  =  T/10  (cm) 

(^c^ln  ) -r   ,  inelastic  scattering  (core) 
C*%n))t   ,  inelastic  scattering  (reflector) 
£    'iteration  criteria  for  3%   accuracy 

D,  'diffusion  coefficient,  core,  Group  I 
D\  dlff.  coef.,  reflector 
f,  f  factor,  core 
f *,  f  factor,  reflector 
(f^s)   for  core 
l\z3)   '  for  reflector 
zla  9     for  core  absorption 
Z.J  *>     absorption  for  reflector 
l/uT,  reciprocal  of  Group  I  lethargy  width 


(See  program)   Same  quantities  as  in  cl, 
(See  Program)  , 

(See  program)  , 


II 

H 

« 

♦1 

II 

II 

argy  width- 

Group 

II 

G-roup 

III 

Group 

T    ( them 

tal) 
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CHAPTSR  II 

Numerical  Integration  Fiethod 

.    Kennedy 

i.  Derivation  of  Equations* 

The  following  diffusion  equation  is  assumed  to 

apply: 

1      o  i_  ^    ^ 

T? V  '  c^  -  5"  .   tf  =  -  5?^  H  Pi   (i) 

■}  2.tr±  ^i   ^  absi  >*i    = = {±) 

d  loge  e 

The  subscript  i  refers  to  a  particular  energy  group.  In 

this  case  (considering;  six  energy  groups)  i  ranges  from 
1  to  6,  where  group  6  neutrons  are  of  highest  energy, 
group  1  of  lowest  (thermal). 

For  grouo  6  equation  (1)  is  modified  to  include 
a  fission  neutron  source.  All  fission  neutrons  are 
assumed  to  be  born  in  group  6. 

equation  (1)  i3  to  be  made  dimensionless,  con- 
verted to  a  standard  form  (to  allov  the  computer  to  use 
one  standard  calculation  form  for  all  groups)  and  integrated. 
It  is  first  multiplied  by  (Z.  ref/IEref )  ,  where  jETref  is  a 
reference  macroscopic  cross-section,  here  taken  as  the 
thermal  transport  cross-section  of  the  reflector. 
Equation  (1)  becomes: 

This  derivation  is  based  on  a  method  developed  by  a.  S. 
Thompson  (J.  Appl.  Phys .  22,  1223)  with  modifications  intro- 
duced to  include  leakage  of  the  fast  group  (6)    and  in  order 
to  exploit  comnuter  advantages.  The  original  method  was 
^veloped  for  hand  calculations. 


-36- 

1  Iref        2  labs  Iref  .  J*i 

3lref        Itr.7     h'  Iref         ^     *  -»W  geT^^) 

tf  logeE 

For  convenience  the  following  diraenaionless  cross-dections 
are  defined: 

ai   s    Zref/Str^ 

bj^  =    Znb3±/  £ref 

di  =    X*±/Xret 

et  ■  Jfi/  Iref 
With  the.se  definitions  equation  (2)  becomes: 

*i— V20  -b  J    0  -  -^fdi^ref  H         (3) 

Zref  ^    x  rei'  *       £  log  L 


__a- 
3 


Next,  to  simplify  the  energy  variable  the  following  d 
nit ions  are  established: 


£>   a  loge  (£fi3s/^i) 
lo"e  (EfiasAth) 

c  =  (loge  Efi3SAth)" 
Thus, 


and  the  operator  (-c  0/  £)  $  )  nay  be  substituted  for 

the  onerator  ( "^ /<&  log  E)  in    ition  i3)>  giving: 

e 
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^Idic£ref01 

ref 

2 


3Z_f  V    *i        t>iZref'°l  a? 


Finally,   the   operator  V    »  which  la  given  in  spherical  co- 
ordinates   (with  radial  variations   only)    by 


sZ 


v2z  a/ar2-  +  f% 


>r 


becomes,  upon  considering  a  dimensionless  spaoe  parameter 
r/R  (where  R  is  the  total  radius,  including  core,  reflector, 
and  extrapolation  distance): 


or  o  -s 

Thus    Equation    (k)    becomes 

3  2ref  R  ^  $ 


(5) 


2      2 

or,    since  the  operator   (R  V    )    in  r-epace  la  equivalent   to 

the  operator   (V     )   in  r/R-apace: 

3ZrefR2  "a?  (6) 

Upon  multiplying  Equation  (6)  through  by  3(  2     f R  ) »  and 
defining 

Al3(  2refR)2 


Equation  (6)  becomes: 

11    1   X     b§ 
Note  that  the  parameter  A  depends  on  the  total  radius  of 

the  reactor  and  is  thus  not  originally  known.   Instead  it  is 

as3umed  (based  on  one— group  calculations)  and  later  corrected 
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if  the  assumed  radius  is  sab-  or  super-critical. 

The  derivative  on  the  ri^ht  hand  3ide  of  Equation 
(7)   way  bf       L   ted: 

so  that:    r±iy2J±   -  b± A  ^±  -  *£C<  Vi        =  -  *£C<V/  <k;   ( ■■) 

The  left  hand  side  of  equation  (8)  represents  loss  from  the 
ifth  energy  group  by  1<    ;e,  absorption,  and  slowing  down, 
ritcht  hand  side  represents  the  input  to  the  i'th  e;roup 
consisting  of  neutrons  slowed  from  the  next  higher  «;roup. 
or  i=l  (thermal  2;roup)  Equation  (£)  is  Modified 

to:       aiY^i  "  biAy*l  =  ^  cd^2  (8a) 

since  no  neutrons  arc  lo3t  from  the  thermal  group  by  slowing 

down. 

The  input  term  in  Equation  (8)  is  modified  for 

energy  group  6  (fast)  to  include  fission: 

*6?\  -   b6^6  "  *  *cA^(>  '     AVei°^  +  £*Vl  iA  (9) 
(The  fission  of  thermal  neutrons  must  be  weighted  by  an 

energy  factor  since  all  other  fluxes  in  (9)  :re  per  unit 
energy  interval  and  must  be  multiplied  by  this  en« rgy  inter- 
val to  become  the' total  flux,  rhe  equation  is  simplified  by 
dividing  tl     l  by  the  energy  interval  with  the  result  that 
it  appears  as  a  denominator  in  the  thi     term.; 

aquations  (8),  (ga; ,  and  (9J  are  of  the  form 
(a^2  -  g2)^  r  h 
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where  g  combines  all  lo3s  te.rr.13  except  leakage,  and  h  repre- 
sents an  input  term,  ror  example,  considering  energy  group  1, 

g2  is  p;iven  by  lb  ?0  and  n  DY  (-Af*^^)  *T°m  (^ai« 

If  this  form  of  Equations  {&) ,  (8a) ,  and'  (9)  is 

L  2 

multiplied  through  by  a^{r/d)      and  the  operator  a.V 

o 

replaced  by  its  equivalent  (a.j   1    d/d(r/AJ  (r/ii)   d/d(r/;j 

( r/n. ) 2 
the  following  expression  (vith  x  -  r/ti)    results: 

aiX4  «T(l/x2)d/dx(aix2)d/dx  -  g2 J  y^  =  ajxMi 
or,  with  ui  =  aj_x2  d/dx, 

(D±2  -  aix4r;2)  ^i  =  atx4h  (10) 

Lquation  (U)  is  of  the  type 

(D2  -  ax4g2)  Y  =  Z 
vrhich  may  be  factored  into 

(D  +  f)(D  -  f )  T  «  Z  (11) 

provided  f  satisfies 

Uf  +  f2  -  ax4g2  =0  (12) 

This  arises  from 

(D  +  f)(D  -  f )  Y  =  Z  =  (D2  -  ax4g2)  Y 

D2Y  +  fDY  -  D(fY;  -  f2l    =  D2Y  -  axV,2Y 

D2Y  +  fDY  -  (Df)Y  -  f(DY)  -  f2Y  *  U2Y  -  ax/fq;2Y 
or,       (Df)Y  +  (f2)Y  -  (ax4g2)Y  =  0 
which  ^ives  Equation  (12)  specifying  f. 

In  Equation  (11)  let  (D  -  f')Y  a  U,  so  that 

(9  +  f)U  =  Z  (13) 
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where     (D  -  f)Y  2  U  (l*+) 

Integration  of  Equation  ilk)   gives 

Y  Z  expifo**   dx«/ax»2)  I    £*exp(-j£*t   dx"/ax"2)U  ■  d*'"  +  (^ 


ax' 

(15) 


and  from  Equation  (13) 


which  combine  to  give  (writing  m  -  /   f  dx'/ax1  ) 

Y  =  m   /*(l/m2)  [/^'"m  Z  dx»  •/**' ,2  +  C21  dx« •  »/ax«  ' ,2  ♦  C 


1 
(17) 


The  constants  C]_  and  C£  are  determined  by  the  conditions 
that 

yl  ;  0  at  x  -  0  and 
Y  =  0  at  x  =  1 
which  is  equivalent  to  requiring  neutron  flux  symmetry  at 
the  reactor  center,  and  requiring  the  flux  to  vanish  at 
the  extrapolated  radius. 

By  the  first  condition  Oz   must  be  zero.  By  the 
second  0]_  is  given  by 


dx"Ax'»2 


01  Z  ~  fo1    (1/m2)  I/*"*  Z  dx'/ax'1 

so  that  Equation  (17)  becomes 

Y=  mf^x[(l/m2)  fx'n  Z   dx»/ax|2  ]dx"/ax««2j         (18) 

Solution  of  Equation  (18)  requires  prior  calcu- 
lation of  f  from  Equation  (12)  with  the  Initial  condition 
that  f(0)  =  0  from  (ll)  and  (10).   Equation  (12)  is  written 
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in  terns  of  differences: 

ax2d/dx(f)  ♦  (f)2  -  ax^g2  =  0 
or,  with  subscript  1  indicatin     ce  point  x  and  subscript 
2  indicating  space  point  (x  ♦  Ax): 

(f2  -  fi)/Ax  =  2(^22X22  ♦  Sixx2)  -  i(f22/a2x2'-  4  f^/a-^2) 
iolving  by  quadratic  formula  gives 

f2  =  -a2x2'^/A  x  +  a2*x2\  a.rx22/^X)2  +  2^l/Ax  ♦  ^2X22  + 

*  *ixi2  -  fi2/aixi2  }  *       (19) 

tors  in  Equations  il#)  and  (10)  are  summarized  for 
clarity: 

g2(loss) 


Energy  Group 

h  ( input ) 

1  (th; 

(-  A  %   Cd2)y*2 

V 

2 

b2A 

b27L  ^  (cd?/A?  )A^ 


(- A^i  cd,.)^  b3A    +  (cd3/Ap    )^S 

(-  A^  cd.)0c  b,X    4  (cd./^e   )^ 

(-  A^  cd6)06  b57V    +    lcd5/Ap     )^ 


V  el<#l  -/Ve^    b6A  4  (cd6/A^  )A$ 
Aft      £ 


i=l 
II.  Fro-     ing  of  Problem  for  Computer 

Whirlwind  I  fWWI)  is  a  digital  computer  capable  of 

performing  simple  mathematical  operations  (addition,  sub- 


tractlon,  multiplication ,  division)  with  extreme  rapidity 
and  under  the  control  of  a  set  of  instructional  called  the 
program,  which  may  be  stored  in  the  computer,  WWI  normally 
operates  with  numbers  that  are  all  less  than  unity,  but,  by 
accepting  a  loss  in  computer  speed  and  storage  capacity,  may 
be  operated  with  unsealed  numbers  (thereby  decreasing 
enormously  the  difficulty  of  setting  up  the  program)  in  the 
Comprehensive  System  (CS  system),  as  is  done  in  this  problem. 

In  the  CS  mode  there  are  1360  registers  available 
for  storage  of  either  instructions  or  numbers.  An  in- 
struction occupies  one  register,  but  a  number  occupies  two. 
Of  the  two  registers  containing  a  number,  the  first  carries 
a  part  of  the  significant  figures  which  represent  the 
number  expressed  in  exponential  form,  while  the  second  car- 
ries the  exponent  and  remaining  significant  figures.  For 
example ,  the  number  123*^  may  be  expressed  as  .123^  x  10-> 
with  .123  stored  in  the  first  register,  and  the  exponent  3 
and  the  remaining  significant  .000*4-  stored  in  the  second. 
Actually,  WWI  operates  with  binary  numbers,  but  an  under- 
standing of  this  feature  is  not  essential  to  the  under- 
standing of  the  problem  since  the  computer  accepts  decimal 
numbers  from  the  programmer  and  returns  decimal  number 
answers. 

The  I36O  registers  which  comprise  the  internal 
storage  capacity  of  WWI  in  the  CS  mode  are  called  the 
magnetic  core  memory  (MCM).   External  storage  is  available 
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in  the  form  of  magnetic  drums  and  tape  but  acquisition  time 
is  large  in  comparison  with  the  time  required  for  the  com- 
puter to  obtain  a  number  stored  in  K«CM  and  is  therefore 
not  used  in  this  problem. 

Any  problem  then  is  limited  to  a  set  of  numbers 

and  instructions  which  occupy  not  more  than  136u  registers. 

In  this  problem  these  registers  are  used  (starting  at  the 

first  register  available,  which  is  number  32)  as  follows: 

Registers  32-20°:  working  storage  (for  numbers 

which  are  calculated  during 
the  problem  and  must  be  tem- 
porarily stored) 
Registers  210-46°:  storage  for  20  value3  of  the 

neutron  flux  in  each  of  the 
six  energy  groups 
Registers  470-509:  storage  for  20  space  points 

(r/R  =  »°5,  .10,  .15  etc.) 
Registers  510-597:  cros3-section  storage 
Registers  5°<--639:  storage  of  constants  other 

than  cross-sections 
Registers  640-end:  storage  for  the  program 

instructions 
Actually  the  program  for  this  ->  roblem  uses  only  some  1^50 
registers  so  that  110  are  available  for  later  expansion  of 
the  problem  if  desire^ . 

For  the  purposes  of  this  problem  there  are  sixteen 


coramands  which  WW I  obeys  in  the  CS  mode.   They  each  affect 
an  operation,  designated  by  three  letters,  and  an  address. 
The  address  may  consist  simply  of  the  number  of  a  register 
which  contains,  for  example,  the  number  which  the  computer 
is  to  obtain  from  MCM.  However,  since  a  programmer  cannot 
know  in  advance  exactly  how  many  registers  to  set  aside, 
it  is  frequently  difficult  to  specify  a  register  number  as 
an  address.   To  eliminate  this  oroblem,  WWI  is  designed  to 
allow  use  of  indefinite  addresses,  called  flads  (for 
floating  addresses)  and  represented  by  a  letter-number  com- 
bination such  a3  al ,  a2 ,  bl ,  b2 ,  etc.   Thus,  if  a  number  is 
to  be  stored  at  the  end  of  the  program  it  is  tagged  with  a 
flad  and  the  flad  used  in  place  of  the  actual  register 
number.   WWI  simply  counts  the  registers  as  the  program  is 
read  in  and  remembers  the  actual  address  of  the  flad.   Then, 
whenever  this  flad  is  referred  to  later  WWI  substitutes  the 
proper  address. 

The  three-letter  portion  of  the  WWI  commands  gener- 
ally calls  for  a  certain  arithmetic  operation.   All  of  these 
occur  in  the  Multiple  Register  Accumulator  (MRA).  For  ex- 
ample, if  a  number  in  storage  is  to  be  multiplied  by  another 
stored  number,  the  first  Is  taken  from  MCM  to  the  MRA.  The 
second  is  then  obtained,  the  required  product  formed  in  the 
MRA  and  then  transferred  to  storage*  Various  operations,  as 
explained  below,  may  either  preserve  or  destroy  the  originally 
stored  quantities. 
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The  sixteen  instructions  used  in  this  problem  are 


the  following: 
Instruction 
Its  al 


lex  al 


ioa  al 


lcs  al 


iad  al 


isu  al 


imr  al 


idv  al 


Meaning 

transfer  the  con- 
tents of  MRA  to 
register  al 


exchange  the  con- 
tents of  MRA  and 
register  al 

clear  MRA  and  add 
to  it  the  contents 
of  register  al 


clear  MRA  and  sub- 
tract from  it  the 
contents  of  al 

add  the  contents 
of  al  to  MRA 


subtract  the  con- 
tents of  al  from 
MRA 


multiply  MRA  by  al 
and  round  off 


divide  MRA  by 
contents  of  al 


Ope.rat.lon 

the  number  in  the  MRA  i8 
transferred  to  register 
al;  the  former  contents  of 
al  are  lost;  the  contents 
of  MRA  are  unchanged 

the  contents  of  the  MRA 
become  the  new  contents  of 
register  al ,  and  vice  versa 

the  MRA  is  cleared  (and 
its  contents  lost)  and  lt3 
new  content  is  the  number 
stored  in  al;  al  remains 

unchanged 

as  above  except  that  MRA 
acquires  the  negative  of 
al;  al  is  unchanged 

the  contents  of  al ,  plus 
the  original  contents  of 
MRA  are  added  and  the 
result  stored  in  MRA;  al 
is  unchanged 

the  contents  of  al  are 
subtracted  from  MRA  and 
the  result  stored  in  MRA; 
al  is  unchanged 

the  number  in  MRA  is  mul- 
tiplied by  the  number  in 
al;  the  result  rounded 
off  to  the  maximum  number 
of  significant  figures 
which  can  be  stored  (about 
8  decimally)  is  retained 
in  MRA}  al  is  unchanged 

MRA  is  divided  by  con- 
tents of  al ,  with,  the 
result  stored  in  MRA;  al 
is  unchanged 
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The  following  two  Instructions  provide  transfer  of 
control  to  the  Instruction  stored  in  the  address  of  the 
instruction.   WWI  normally  operates  sequentially  from  one 
register  to  the  next.   These  two  instructions  circumvent 
this  sequential  operation.   As  an  example,  the  last  in- 
struction in  an  iterative  problem  might  contain  a  transfer 
of  control  back  to  the  first  in  order  to  perform  the 
iteration. 
Instruction      Meaning  Operation 

isp  al        control  transfer   take  the  next  instruction 

from  register  al  and  pro- 
ceed from  there 

lcp  al       conditional       as  above,  if.  the  contents 

control  transfer   of  MRA  are  negative;  other- 
wise, ignore,  and  take  the 
next  Instruction  In  sequence 

The  next  five  instructions  have  to  do  with  operations 
performed  with  a  counter.   In  WWI  a  counter  may  be  set  aside 
to  allow  repetition  of  an  instruction,  or  set  of  instructions, 
a  specified  number  of  times.   The  specified  number  is  the 
criterion,  the  number  of  times  the  operation  has  been  per- 
formed is  the  index.   WWI  sets  aside  both  criterion  and 
index  (with  index  Initially  zero),  cycles  the  index  once 
each  time  the  operation  is  repeated  until,  upon  comparison 
the  index  is  found  to  equal  the  criterion  (or  if  the  index 
exceeds  the  criterion,  as  pointed  out  below)  at  which  time  the 
computer  proceeds  to  the  next  instruction  in  sequence. 
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Instructlon       Meaning:  Operation 

lso  J         select  counter  J   select  counter  J  for  fol- 
lowing operations  with 
counter 

(selecting  a  counter  J,  by  lsc  Jf  does  not  In  Itself  set 

criterion  and  index;  this  must  be  done  by  an  additional 

step  prior  to  the  first  use  of  counter  J.) 

Instruction       Mean  inn:  Operation 

icr  a        cycle  reset        set  criterion  equal  to  m, 

index  zero 

lot  al        cycle  count_        increase  the  index  register 

by  one  and  compare  with 
criterion;  if  the  index  is 
less  return  to  register  al; 
otherwise  take  the  next  in- 
struction in  sequence 

ici  n         cycle  increase     Increase  the  index  by  n 

(this  may  3et  index  greater 
than  criterion,  but  is  use- 
ful in  operating  with  more 
than  one  counter) 

led  n        .cycle  decrease     decrease  index  by  n 

(Note  that  the  counting  instruction,  let  al ,  is  the  only 

counter  ^™..mand  that  transfers  control.   It  is  therefore 

usually  the  last  instruction  in  an  operation  which  is  to  be 

repeated. ) 

During  operation  with  a  counter  WWI  may  treat 

numbers  (or,  in  special  cases,  instructions)  which  are  stored 

in  sequence  if  +  c  is  appended  to  the  address  of  the  first 

number  to  be  treated.   For  example,  to  place  ten  numbers 

whioh  are  stored  consecutively  in  ten  other  consecutive 

registers  it  suffices  to  select  a  counter  with  criterion  ten, 
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iadex  zero,  obtain  and  transfer  the  first  number,  with  +c 
appended  to  the  address  section  of  the  obtain  and  transfer 
instructions,  and  transfer  control  back  to  the  first  instruc- 
tion (in,  say,  register  bj_)  with  an  ict  bl.  Lach  time  the 
computer  notes  the  +c  in  the  address  it  set3  c  equal  to 
twice  the  index  of  the  counter  in  use  and  thus  transfers 
the  ten  desired  numbers.  The  computer  set3  c  to  twice  the 
index  because  of  the  fact  that  the  numbers  occupy  two 
registers.  Phis  however  is  not  the  case  if  a  +c  is  added 
to  the  address  section  of  an  iso  command.  Here  the  address 
mujt  be  -mother  instruction,  occupying  one  register,  so 
the  computer  conveniently  sets  c  equal  to  the  inde: , 

The  final  instruction  is  used  in  connection  with 
subroutines  within  the  program.  Identical  operations  which 
are  to  bt  performed  at  different  points  in  the  program  nay 
be  set  an  as  a  subroutine  in  order  to  avoid  duplication  of 
instruction  storage  each  time  the  required  calculation  is 

i.  Then,  at  each  point  the  d     itine  is  needed,  an  isp 
instruction  transfers  control  to  the  subroutine.  However, 
to  leave  the  subroutine  and  return  to  the  proper  ooint  in 
the  main  program  an  isp  instruction  at  the  end  of  the  sub- 
routine must  have  the  address  of  the  register     iiately 
following  the  register  which  originally  transferred  control 
to  the  subroutine.  This  is  accomplished  bv  the  final  instruc- 
tion to  be  considered  here. 
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Instractlon       Meaning  Operation 

ita  al        transfer  address   the  address  section  of  the 

Instruction  in  al  is  ohanged 
to  the  register  immediately 
after  the  register  con- 
taining the  latest  lap  in- 
struction which  has  been 
obeyed 

For  example,  to  obtain  the  square  root  of  a  number  at  several 
points  in  the  main  program,  a  square  root  subroutine  may  be 
set  up  at  the  end  of  the  program.   Then,  each  time  the 
operation  is  to  be  performed  the  number  is  placed  in  the 
MBA,  control  transferred  to  the  subroutine,  the  root  calcu- 
lated in  the  subroutine  and  placed  In  the  MRA,  and  control 
transferred  back  to  the  main  program  to  the  register  im- 
mediately following  the  original  transfer  point. 

In  the  subroutine  the  first  instruction  would  be 
ita  bl,  with  bl  the  last  register  in  the  subroutine.   In  bl 
would  be  stored  lsp  — •   Thus  the  ita  instruction  would 
cause  transfer  of  control  back  to  the  main  program  at  the 
proper  point. 
Ill*   Programming  of  Specific  Problem. 

The  problem  to  be  solved  here,  derived  in  Part  I 
of  this  chapter,  may  be  stated  as 


0  ax1 


2 


Zda 


hi 


ax1"2 


where  x  ia  a  dlmenaionless  apace  parameter  r/R  total 
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a  ij  a  stored  dimensionless  cros3-section; 
Z  is  the-  source  term,  representing  fission  neutrons 
(for  the  highest  energy  group)  or  slovred  neutron3 
(for  all  lover  energy  groups);  and 
f  is  calculated  from  a  finite  difference     roxirnation, 
Fundamentally  the  computer  proceeds  by  the  fol- 
lowing    3: 

(1)  Calculate  f (xj ; 

(2)  divide  f(x)  by  ax2  and  integrate  from  u  to 


x  bv  trapezoidal  summation; 

(3^  Calculate  and  store  exp  / 

yo 
(4)  Calculate  source  term  Z,  initially  from 


an  assumed  flux  and  critical  radius; 

(5)  Divide  Z  by  ax  and  multiply  by  the  expo- 
nential computed  in  step  (3),  preserving  the 
stored  exponential  for  later  r>  -use; 

(6)  Integrate  /   |(Z/ax  )(exp  /  f  dx/ax2)  dx  ; 

(7)  Divide  the  above  integral  by  the  square  of 

the  exponential  (still  preserving  the  expo- 

2 
nential)  and  by  ax  ; 

(V)    Reverse  the  order  of  storage  of  the  quantity 

determined  in  step  (7)  so  that  the  integral 

from  1  to  x  may  be  determined  (this  step 

is  necessary  since  v.M  counters  increase 
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rather  than  decrease  following  each  repetition). 

Step  (8)  effectively  cancels  the  minus  sign  which 
appears  in  the  input  term  to  each  group,  with  the  result 
that  the  calculated  neutron  flux  is  positive,  as  it  should  be, 

(9)  Reverse  the  order  of  storage  of  the  integrals 
determined  in  step  (8)  and  multiply  by  the  exponential  term, 

(10)  Store  the  results  of  step  (9)  and  proceed  to 
the  next  energy  group  calculations ,  changing  f,  a,  and  Z  in 
steps  (l)  through  (9)« 

Following  the  sixth  application  of  these  steps  the 
new  thermal  flux  has  been  computed;  WWI  proceeds  then  to  the 
next  step: 

(11)  Compute  a  new  critical  radius  parameter  ^': 


A': 


A 


where  the  primed  quantities  are  the  newly  computed  values. 

Step  (11)  is  based  on  the  fact  that  once  the 
proper  solution  for  the  flux  is  determined  0.,  will  not 
change  with  more  iterations,  as  long  as  the  radius  is  equal 
to  Rcrit*   ^n  ^e  other  hand,  if  the  radius  assumed  was  sub- 
or  super-critical  0+^   will  be  less  than  or  more  than  0^ 
(assumed)  respectively. 

(12)   Compare  the  two  values  of  the  radius  parameter; 
if  they  differ  by  more  than  an  acceptable  error  recalculate , 
based  on  the  new  values  of  flux  and  radius;  if  the  reactor 
Is  critical,  as  evidenced  by  a  specified  small  variation  in  7^ 
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print  out  specified  information  (e#g«,  flux,  radius,  power 
distribution)  and  stop. 

The  comouter  may  be  stopped  at  the  end  of  a 
problem  by  an  13T0P  instruction.  Alternatively,  if  more 
problems  are  to  be  considered  they  may  be  listed  after 
the  firs-     hlpm  and  the  computer  stopped  with  the  special 

■t ruction  sp26.  when  it  reaches  this  command  at  the  end 
of  the  first  problem  the  comeuter  automatically  reads  in 
the  remaining  problems  consecutive] v.  in  this  way  a 
considerable  savins;  in  time  may  be  realized,  since  manual 
effort  \  otherwise  required  to  start  each  problem]  i3  mini- 
mized and  since  many  of  the  constants  which  arc  already  in 
storage  will  not  change  from  problem  to  problem,  thus 
decreasing;  the  time  required  for  the  computer  to  convert 
the  input  decimal  numbers  into  binary  numbers.  After  the 
final  rroblera  an  ioTUF  command  may  be  read  in  to  stop  the 
computer.  Otherwise  the  operator  will  stop  the  computer 

Hially  when,  upon  completion  of  the  last  problem,  the 
computer,  under  the  sp26  instruction,  calls  for  more  infor- 
mation. 

as  mentioned  in  part  II  of  this  chapter,  the 
numerical  values  of  this  problem  are  storer  in  the  first 
portion  of  MCM.  The  first  instruction  is  stored  in  register 
ent  of  the  program  for  this  problem      start 
at  this  point  with  numbers  specified  as  they  appear. 


Instruction 

alf IscO 

icrppl 

iscl 

icrpp2 

isc2 

icr20 

isc3 
icr20 
isc4 
icr6 
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Operation 

Sets  counters:  0  Tor  iteration  over 
core  region,  3  over  reflector,  2  and 
3  over  entire  reactor  (20  space  points) 
and  4  over  six  energy  groups;  ppl  and 
op2  are  preset  parameters  which  are 
set  at  the  op^  rung  of  the  program  in 
orde-  that  the  rest  of  the  program 
may  remain  unchanged 


ica636 
its598 

isc2 

a3,ica430+c 

imr638 

its430+c 

icta3 

f2,isc0 

a2,ica510 

imr254+c 

its3^+c 

ica512 

imr298+c 

iad3^+c 

its34+c 

ica51^ 

imr342+c 

iad34+c 

itsS^+c 

ica5l3 

imr386+c 

iad34+c 

its3^-K: 

ica5l8 

imr430+c 

iad34+c 

imr598 

imr470+c 

imr^70+c 

its3^+c 

icta2 


Stores  assumed  radius  parameter  in  598; 
636  contains  latest  calculated  A- 


Scales  0 

value 

register 


to  any  desired  center-line 
(value  initially  set  in 
63S) 


th 


Select's  counter  0  to  iterate  over  core; 
double-length  registers  510  through  51S 
contain  fission  parameters;  254  through 
46B  contain  fluxes ;  input  term  Z  to 
fast  energy  group  i3  determined  as  the 


sum  of  each 


[l>Jtf  $i ) 


(A)  in  59^,  x^  (x  in  470  through  508) 
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lca574 
Its6l8 

ica57o 
it3620 

Ica530 
its622 
Ica532 
its624 


J tores  loss  parameters  (b*]|  cd/t$  )   and 
a  (^rcf/?trj  of  group  6  in  workin  ;  r 
at era  618-6^4  (this  is  done  in  order 
that  other  group  parameters  may  later 
be  placed  in  th>-  ;   same  working  regi- 
sters arid  the  remainder  of  the  program 
used  for  all  energy      sj 


ispml 

a4,Iac4 
ictml 

ispx6 

ml,iscO 

a5,ica6l8 

imr470+c 

imr470+c 

its78+c 

icta5 

iscl 

,ica620 
imr470+pp3+c 
imr470+pp3+c 
lts7o+PP3+c 
icta6 

ica620 

iad6i8 

imr6l4 

imr470+pp3 

Imr470+pp3 

its78+pp3 

iscU 

a7,ica622 

imr470+c 

imr470+c 

imr598 

itsl22+c 

icta7 

iscl 

,ica624 
imr470+pp3+c 
imr470+pp3+c 
imr598 

Itsl22+pp3+c 
icta8 


causes  a  jump  around  the  next  3  instruc- 
tions which  provide  entry  for  groups  1-5 
a4  is  flad  to  label  entry  point  for  groups 
1-5;  next  register  jumped  unless  new  $.. 
has  been  calculated 

jumps  to  0£n  subroutine  in  x6  when  new 
$tn  has  been  computed 

2 
calculates  ax  for  cor(  region 


as  above,  for  reflector  region;  pp3  is 
set  at  the  opening  of  the  program  to 

Lee  (because  oi   double-lengt]     i- 
StersJ  the  estimated  value  of  Rcore/R..  ; 
thus,  in  order  to  specify  a  reflector 
core  combination  it  is  only  necessary 
to  change  ppl,2,3i  not  program 

places  av<  ragi  value  of  ax  at  cor  - 

Elector  interface;  614  contains  0.5» 
used  i    ently  in  averaging 


calculates  loss  tern  (  te.  .  ,  b^.+^jjcd/k^  ) 
times  x2  for  core 


as  above |  for  reflector 


lea 

>4 

imr470+pp3 

irai»598 
Ital2&fpp3 

ica626 

idv47 
idv470 

Imr 

l3pE2 

its 
ica 
idv470 

lad 
itsl68 

isc5 

icr 

imr47C 
idv78+c 

lad 
lad 

■■ 
iad3i  i-fc 

idv47G 
idv47 

iad:l22M-c 

iad!22+c 

imr8o+e 

ispz2 

itsl70+c 

icsoO-fc 

idv^TO 

lad ] 70+c 

its:7r+c 

icta9 

isc 

r   ,  Lcal6 
Idv78+c 
its i 60+c 
ictgl 
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as  above,  at  core -re flee tor  interface 


places  zero  in  ,   20 

Initially  set  to  zero  and  le 

to  .. 


icui 

"is ion  by  zero) 


zn  is  souaro  rout  subroutine 


calculates  remaining  19  values  of  f 


calculates  f/-3.xz 


isc2 

bl,ical'56+c 

iadl68+c 

Imr6l4 

imr^70 

it3l22+c 

ictbl 

isc2 

b2,ical22+c 

iadl20+c 

itsl22+c 

ictb2 

isc2 

bjj,ical22+c 

ispzl 

itsl22+c 

Iotb3 

Isc2 

b5#ica3^+c 

imrl22+c 

it334+c 

ictb5 


-5$- 

forms  areas  of  f  dx/ax*  prior  to 
trapezoidal  integration  (/-/Aj  sAx) 


sums  above  areas  in  trapezoidal  inte- 
gration P20  prpviouslv  set  to  zero) 


calculates  exp( 
zl  identifies 

onential  subroutine 


fx. 

I      f  dx/ax  )  s 
Jo 


flad 


source  term  Z  times  exponential 


isc2 
bo,lca32+c 

iad3^+c 

imr6l4 

imr^O 

itsl68+c 

ictb6 

ica626 

it3l56 


forms  areas  of  (Z  times  exponential) 
p^ior  to  trapezoidal  integration 


iets  166  to  zero 


Isc2 

b7,ical68+c 

iadl66+c 

itsi68+c 

ictb7 

iscO 

fi,ical68+c 

idvl22+c 

idvl22+c 

i'dv^70+c 

idv470+c 

idv6l8 

itsl68+c 

ictfl 


sums  areas,  in  trapezoidal  integration, 
to  compute  /x 

/      Z    (exp)    dx/ax 
Jo 


for  core    region: 
rx 

Z    (exp)    dx/ax*1 
'o 


f 


( exp )      ax 


Ical68+pp3 
its74 

iscl 

g3,ical68+pp3-»-c 

idvl22+pp3+c 

idvi22+pp3+c 

idv470+pp3+c 

idv470+pp3+c 

idv620 

itsl68+pp3+c 

ictg3 

Ica620 

Iad6l8 

imr6l4 

iex74 

idvl22+pp3 

idvl22+pp3 
idv470+pp3 

idv470+pp3 
Idv74  , 
Itsl68+pp3 

isc3 
icr20 

leU9 

kl,Isc3 

ical68+c 

icdl 

Isc2 

its34+c 

ictkl 

Ica626 
its74 

isc2 

. ica36+c 
iad34+c 
imr6l4 
imr470 
its7o+c 
ictbo 

ica626 
its7o 
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ts  aside  (in  74) 
at  core-reflector 

for  reflector  region: 
/x 


(exp)  dx/ax 
interface 


2 


Z  (exp)  dx/ax' 


(exp)   ax^ 


as  above,  at  core-reflector  interface 
(usin?  n     tor  set  aside  in  11) 


reverses  storage  of  above  (over  entire 

ctarj  in  order  to  integrate  from 
1  to  x 


seta  register  ?/  to  zero 


forras  areas  prior  to  integration 


sets  register  76  to  zero 
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isc2 

b9,ica78+c 

iad76+c 

its78+c 

ictb9 

isc3 
icr20 

Icil9 

k2,isc3 

ica76+c 

icdl 

isc2 

its3^+c 

ictk2 

isc2 

k3,ica3^+c 

Its7B+c 

ictk3 

isc2 

cl,ica7#+c 

Imrl22+c 

its34+c 

ictcl 

isc4 

isppl+c 

pl,ispxl 

Ispx2 

ispx3 

ispx4 

ispx5 

ispa4 

x6,isc2 

Ica626 

e2,iad34+c 

Ictc2 
itsll6 

Isc2 

ica626 

c4,iad430+c 

ictc4 

idvll6 

isu508 

imr 

iad508 
imr598 
its636 


sums  above  areas  in  trapezoidal  inte- 
gration 


returns  above  integrals  to  original 
order  of  storage 


(exp)  dx/ax2 


(exp) 


dx 


ax 


counter  L,   selects  energy  group  by  jump 
to  xl  through  x5  (which  store  0^  through 
02  and  position  constants  for  next  co? 
puting  tyr}   through  0i);  the  isp  a4 
instruction  resets  counter  L   and  pro- 
ceeds to  the  0i  storage  routine 


to  determine  radius  criticality  calcu- 
t.es  an  approximation  of  Z1  0L  tx)  dx 

based  on  the  newly  computed  $i   in  the 
registers  following  34 

as  above,  based  on  previously  computed 
(or  assumed)  0^ 

computes  ratio  of  above  integrals  if 

computes  f -   1 

multiplies  (tf-  1)  by  any  desired  factor 

control     iity  of  cor     .on,  plus  1 
computes  new  ^-  :    r 

/t,  =  \   (  i+\r-  #n) 


isuj 

iad3^ 

its!22 
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lin-  of  ^f. 

at   reactor   center  to   non    ilj 
Lculated  ^t^ 


lsc2 

,Ica3^+c 
idv!22 
lmr638 

ita430+c 

ict 


nur  /'         o  corres« 

pond  sSn    i      (in 

at  center 


store  normalized  ^th  in  assigned 
registers 


idvl 

ispz2 

iM0A+i23.^56c 

isu- 
idv 
iad( 

icpal 
lca! 

lad. 
icpal 


calculate  it-total  from    ~A> 


print  out  Htota]_  (sample  number  given 
to  specify  number  of  significant 
figures  desired.) 

determine  critic ality  by  comparing 
the  per  centage  change  in  A 
with  an  acceptable  error  stored  in 
register  604;  if  error  is  too  large — 
either  plus  or  minus — return  to  al 
(the  original  starting  point  of  the 
program)  and  recalculate;  if  error  is 
acceptable,  continue  to  next  sequence 
of  instructions 


t 

+0 
+20 


ecifies  format  of  solution:  blocks 
of  20  numbers  arranged  5  per  line , 
with  a  tab  between  numbers 


illOA+i  , 
ict 


prints  out  2u  values  of  the  calculated 
flux  for  each  of  6  energy  groups  con- 
sidered 


i3G 

OA+i.23^5f 

ict. 
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z2 
j  L 

ictc 

lsc2 
ciO,ica« 

•i.23k 
ictclO 

1bc2 
eil,ica25^+c 

JA-H.23V 
Ictcil 

isc2 

012,10.1210+0 
1.2  3Ji 


cO 

-520 

J+G 

•i-C 

itc3Jl. 

ica52'4 

lmr3'L- 

iad3ir+c 

its3^+c 

ica^ 

inr2J 

lad3^+c 

its3'++c 

icar 

imr25^+- 
iad3^+c 

its3^+c 

lCt: 

ica3^ 
1bu36 
lad  34 

lts!22 


calculate  the  conribution  to  the 
total  power  distribution  from  each 
group  (fission  cro33-3ections-- 
ighted  by   A^  ,/c  for  ^  3  l   5 — 

stored  in  520,522 ,£24,526,52$) 


linear  extrapolation  to  determine 
power  at  center 
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iscO 

,ica34-»-c 

idvl22 

..J 
Ictm7 


,+C 

lt6210+c 
let    . 

ia578 
tts6i 
lca530 
Its620 

ica53^ 

its62 

ica536 

its62;'i 

iscO 

d6,ica3^+c 

lmr470+c 

imr!l-70+c 

lmr598 

its3JM-c 

ictdo 

ica3^-*-PF 

its602 
iscl 

e2,ica34+pi 
imrVfO+ppj/hc 
imr-V/  re 

ict 

lca55J* 

Iad55 

imr 

Imiv 

irar470+p 

imi470+pi 

its3]++P, 
isp 


print  out  power  distribution 
normalized  to  unity  at  center 


insert  stopping  command  (e.g., 
iSTOP  or  sp26,  as  discussed 
on  page  54; 

subroutine  to  be  used  after 
calculation  of  $,:    stores  0,- 

in  assigned  registers,  places  the 
parameters  of  the  next  lover 
energy  group  in  the  working 
registers  and  calculates  the 
source  term  Z  for  the  next 
lover  group;  reenters  program 
at  entry  point  set  aside  on 
page  56 


(THE  FOLLOWING  PAG:     fTAINS 
FU'u    IMILAR     OUTINES  WHU 

■:tion  aftfr 


CALCULATION  OF  J2L ,   0 


5'   >V 


fyt: 


) 
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x?,isc2 
d2,ica34+c 

:lts254+c 
ictd2 
icn582 
lts6_ 
,ica534 
its620 
ica5. 
itco22 
ica540 
its624 
iscO 

d7,ica34+c 
lrar470+c 
ir470+ 

lmr55 

s34+e 
1g  td? 

ica34+pp3 
its602 
iecl 

e3*ica34+pp3+c 
imr470+op3+c 
imr470+pp3*c 
imr593 
imr56o 
-s34+pp3+c 

icte3 

ica558 
iad56o 
Imr6l4 
imr602 

Imr470+pp3 
imr470+pp3 

its34+pp3 
I  spa  4 

.  i,  exponential 
j  square  root 
i  START  AT  640 


x3,isc:^ 
<33,ica34+c 
its293+c 
ictd3 

ica^c 

>i 

ica588 

0620 
ic 

its622 
ica544 
its624 
iscO 

;ica34+c 
Imr470+c 

-470+c 
lmr598 

ir562 
ito34+c 
ictd8 
ica34+pp3 
its602 
iscl 

«  r,ica34+pp3+c 
Imr470+pp3+c 
iittp4704oo3+c 

Imr59 

imr564 

Its34+pp3+c 

icte4 

ica5o2 

iad564 

imr6l4 

imr602 

Itar470+pp3 

imr470+pp3 

imr5j)8 

Its34+pp3 


;isc2 

,ica34+c 
its342+c 
ictd4 
ica530 
its6l8 
iea., 
its620 
ica5'r 
Its622 
Ica548 
tts624 
iscO 

ij,lca34+c 
imr470+c 
imr470+c 
imr5  .■ 

its34+c 

ictd9 
ica34+pp3 
lts602 
iscl 

Lca34+pp3+c 
imr4?0+pp3 
Imr470+pp3+c 
imr598 
imr568 
its34+pp3+c 
icte5 
ica^oo 
iad568 
Imi*6i4 
imr602 
irar470+pp3 
imr470+pp3 
imr598 
Its34+pp3 
isoa4 


x5,isc2 

d5»ica  34+c 

its386+c 

ictd? 

ica594 

its 

Icafj 

its:. 

ica550 

its. 

;24 

iscO 

ei,ica34+c 

lmr470+c 

imr470+c 

inir! 

imrf  7 

its34+c 

ictel 

ica34+pp3 

its602 

iscl 

i  -j,ica34+pp3+c 

imr47C+pp3+c 

imr470+pp3+c 

imr598 

imr-572 

its34+pp3+c 

icte6 

ica570 

iad572 

imrul4 

irar6Q2 

imr470+pp3 

imr470+pp3 

imr598 

itB3.4+pp3 

a4 


ispa4 

subroutine  (LSR  PUl) 
subroutine  (LSR  FU  2b) 

(terminates  read-in  and  instructs  computer 
to  take  first  instruction  from  register  64O) 
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(Al]     wen-numbered   r  •  >   include,    the   follow ;  id 

nura  provide  double  ;c) 

;er(s)  bored  Da 

32-  Initial         do  as  storage  during 

■  gram  cal<::     "ons  by  insertin 
roes 

430-  Insert  2         >f  assurr        lal 

flux 

Insert  20 
(x*r/R».05,.iO,.15j 25,1. 

tsion  para      :  "v'^^jj.j 

■  4  c 
»  3c 

^■f     lc 

11 
Are/  *A^2>ef 


532 
534 

>,  except  f  »r 
P ) 

554  In 

in 

572 


It 

II 

H 

5r 

II 

4c 

II 

4r 

11 

II 

3r 

II 

2  c 

Z^lr, 

1  ^r 
*  lc 

11 

cM?^- 

ii 

11 

—  _ 

11 

90 

5r 

11 

M 

11 

3o 

u 

2r 
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ansionlea  ,'c 

7H 


!3BC 

ic 

3c 

2c 
1c 

/\-    (ini 
^  A) 

o 

t 
In  A    "    C 

■ 
-   10 
1  .  A 

lb— 6  :rs 

i 

i 
■ .  l 
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FLOW  SHEET  FOR  SIX  GROUPS 
(Superscript  i  denotes  energy  group,  from  1  to  6) 
(Subscript  n  denotes  space  point,  fron  1  to  20) 


^Uiy 


Position  Group  i 
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FLOW   SHEET  FOR   SIX   GROUPS    (Continued) 
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FLOW  SHEET  FOR  SIX  GROUPS  (Continued) 
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CHAPTSR  III 
R53ULT3 

!•   Reactor  Compositions  Studied 

a»   Core  Materials »   In  order  to  simplify  the 
problem,  a  choice  was  made  as  to  the  type  of  fuel  element, 
coolant  and  moderator.   The  choice  was  based  on  a  considered 
adaptability  to  high  output  mobilo  reactors.   High  power  re- 
quirements, fioupled  with  a  need  for  light  weight  and  small 
size,  dictated  the  use  of  highly  enriched  uranium  and  clad- 
ding materials  which  would  withstand  high,  temperatures. 

The  fuel  element  waa  a  dispersal  element  of 

uranium  dioxide  in  a  matrix  of  stainless  steel.   3ome  calou- 

1 

lations  have  been  made  on  this  type  of  fuel  element  ,  and  good 

high  temperature  properties  are  envisioned.   The  amount  of 
uranium  epoxide  in  the  steel  matrix  is  kept  low  in  orde^  to 
minimise  radiation  damage  and  to  improve  heat  transfer 
properties.   A  similar  fuel  element,  with  the  steel  replaced 
by  aluminum,  was  operated  at  the  Genera  conference  on  the 
peaceful  uses  of  atomic  energy. 

The  uranium  considered  is  enriched  to  90$  U-235* 
The  density2  of  U02  is  10,95  gm/cc# 

Stainless  steel  is  also  used  as  the  cladding  for 

3 
the  fuel  element.   Its  choice  is  dictated  by  the  fact   that 

it  will  withstand  sodium  up  to  about  1600  F.   Its  absorption 
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cross  section  presents  little  problem  with  the  fuel  enrich- 
ment considered  here. 

Beryllium  is  chosen  as  moderator  because  of  its 
good  moderating  properties,  light  weight,  and  good  thermal 

properties.  It  is  protected  from  the  sodium  coolant  by  stain- 

3 
less  steel,  since  beryllium  withstands  sodium  only  to  about 

1000  F  and  higher  temperatures  than  this  are  considered. 

Choice  of  sodium  as  the  coolant  is  based  on  its  good 
heat  removal  characteristics  and  high  boiling  point.  A  high 
boiling  point  eliminates  the  need  to  pressurize  the  core. 

The  reflector  is  considered  , to  be  composed  of 
various  proportions  of  Be,  Na,  and  stainless  steel.  Reflector 
cooling  is  necessary  since  some  five  per  cent  of  the  heat 
may  be  liberated  in  the  reflector.  This  amount  of  heat  will 
be  considerable  for  high  output  reactors.  Stainless  steel 
is  required,  as  in  the  core,  to  protect  the  beryllium  from 
the  corrosive  action  of  high  temperature  sodium,  and  for 
structural  purposes. 

With  a  basic  reactor  chosen,  variations  in  the 
proportions  of  the  materials  were  studied  in  order  to 
determine  their  effect  on  the  power  distribution.  Table  III 
of  Appendix  B  summarizes  reactor  compositions  studied. 

The  effect  of  reactor  poisons  was  not  considered. 
All  of  the  reactors  were  assumed  to  be  ho  ogeneous  and  in 
the  steady  state. 
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*>•   Cross  Sections  of  Reactor  Materials •   Tho 
cross  sections  of  the  materials  considered  are  listed  in 
Tables  I  and  II  of  Appendix  B  •   They  are  based  on  data 
from  "Neutron  Cross  Sections/  BNL  325  (August,  1955* 
Geneva  Supplement). 

c«   Reactor  Configuration <>   For  reasons  of  mathe- 
matical simplicity  the  reactor  considered  i3  spherical, 
with  a  reflector  of  uniform  thickness.   Because  of  angular 
symmetry  the  spherical  reactor  may  be  resolved  to  a  one- 
dimensional  problem.   Cylindrical  or  rectangular  geometry 
would,  on  the  other  hand,  require  two-  or  three-dimensional 
analysis.   While  it  may  be  argued  that  spherical  reactors 
will  not  find  extensive  use,  the  results  from  an  analysis  of 
spherical  reactors  can  be  extrapolated  to  cover  reactors  of 
other  configurations. 

The  magnitude  of  the  multi-dimensional  problem 
is  such  that  simplification  to  a  single  dimension  is  made 
frequently.   Other  than  the  spherical  case  considered  here, 
the  only  one-dimensional  geometry  is  that  of  the  infinite 
slab*   Spherical  geometry  Is  superior,  since  it  gives 
results  which  may  be  applied  to  reactors  of  other  configu- 
rations whose  dimensions  along  various  axes  are  not  greatly 
different. 

d»  Fuel  Concentration*  The  proportion  of  fuel  in 
the  reactor  has  an  important  effect  on  the  radius  of  the 
critical  reactor.  It  is  also  related  closely  to  the  power 
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level. 

For  the  mobile  power  reactor  envisioned  here, 
there  are  certain  restrictions  on  the  size  which  must  be 
considered.   It  cannot  be  so  large  that  its  mobility  is 
questionable,  and  yet  it  mu3t  be  large  enough  so  that 
enough  coolant  may  be  pumped  through  to  remove  the  required 
heat. 

From  elementary  heat  transfer  considerations 
(Appendix  A  ) ,  a  core  radius  of  from  two  to  four  feet  was 
found  to  be  sufficient  for  heat  removal,  and  since  this 
size  appeared  reasonable  for  mobile  applications,  the  amount 
of  fuel  was  adjusted  in  order  to  give  a  radius  within  these 
dimensions. 

e«   Moderator  Proportion*   The  amount  of  moderator 
considered  varied  from  zero  to  sixty  per  cent  of  the  core 
volume.   This  had  the  effect  of  changing  the  reactor  from  an 
essentially  thermal  one  to  an  intermediate  one,  in  which  the 
greater  part  of  the  fissions  take  place  at  energies  appre- 
ciably above  thermal. 

f .   Reflector  Composition.   It  wag  determined  that 
as  little  as  two  per  cent  of  the  reflector  volume  was  re- 
quired to  be  3odium  In  order  to  carry  off  heat  generated 
within  the  reflector.   The  amount  of  stainless  steel  was 
arbitrarily  taken  as  five  per   cent. 

By  assigning  the  remainder  of  the  reflector  volume 
to  the  beryllium,  it  was  found  that  the  neutrons  returned  to 
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the  core  were  largely  thermal.   In  order  to  decrease  the 
number  of  moderated  neutrons  reflected  back  to  the  core,  ad- 
ditional sodium  was  Introduced  into  the  reflector  and  the 
beryllium  decreased,  so  that  the  percentage  of  beryllium 
considered  ranged  from  a  low  of  55%   to  a  high  (one  run  only) 
of  95*. 
II  •   Effect  of  Various  Parameters. 

a»  Fuel  Concentration*  Apart  from  its  effect  on 
the  critical  radius  of  a  reactor,  the  fuel  concentration 
plays  an  important  part  in  determining  the  energy  spectrum  of 
neutrons  in  the  reactor.   Neutrons  are  born  fast  and,  during 
the  slowing  down  process,  may  be  absorbed  by  reactor  ma- 
terials. Absorption  during  slowing  down  increases  the  number 
which  must  be  born  fast  in  order  that  enough  remain  after 
slowing  down  to  cause  further  fissions  and  perpetuate  the 
chain  reaction. 

An  examination  of  the  cross  sections  of  the  ma- 
terials considered  here  shows  that  the  predominant  absorber 
of  neutrons  having  energies  greater  than  thermal  is  the 
uranium.   Increasing  the  fuel  thus  Increases  the  flux  In 
higher  energy  groups. 

Actually  there  are  three  things  which  may  happen  to 
a  neutron  while  slowing  down:   it  may  be  absorbed;  it  may  be 
lost  through  leakage;  or  it  may  undergo  scattering  collisions 
and  be  slowed  down.  The  fuel  concentration  Is  so  small 
that  it  can  have  little  effect  on  the  slowing  down  power. 
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The  fuel's  important  effect  la  on  absorption  tendencies  and, 
as  a  consequence,  on  the  leakage,  since  an  increase  in  fast 
flux  with  more  fuel  Increases  leakage. 

It  is  important  to  note  that  of  the  neutrons  ab- 
sorbed by  the  uranium  during  the  slowing  down  process, 
rouerhly  four—fifths  cause  fission  of  the  absorbing  nucleus. 
Thus  more  fuel  has  the  added  effect  of  changing  an  originally 
thermal  reactor  toward  a  reactor  in  which  fissions  at 
energies  above  thermal  predominate. 

Figure  1  demonstrates  the  above  effect  of  fuel  con- 
centration on  the  ratio  of  fast  to  thermal  flux.   The  ratio, 
calculated  at  the  reactor  center,  is  plotted  versus  fuel 
concentration  for  twenty  different  reactors  all  having  at 
least  30$  beryllium  in  the  core. 

It  should  be  pointed  out  that  the  ratio  of  fast  to 
slow  flux  has  an  important  effect  on  the  ability  of  the 
refleotor  to  assist  in  levelling  the  power  distribution.  A 
greater  number  of  fast  neutrons  leaking  into  the  reflector 
increases  the  number  of  slowed  neutrons  which  the  reflector 
can  return  to  the  core  to  compensate  for  slow  leakage,  which 
would  otherwise  depress  the  power  at  the  edge  of  the  core. 

b»   Moderator  Concentration «   As  mentioned  above, 
the  ratio  of  the  fast  to  the  slow  flux  depends  on  leakage, 
absorption,  and  slowing  down.   Changing  the  moderator  pro- 
portion in  the  core  affects  the  first  and  last  terms.   For 
simplicity  of  discussion  its  effect  on  leakage  may  be 
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neglected,  if  the  central  regions  of  the  reactor  are  con- 
sidered. 

In  comparing  the  remaining  two  terms — absorption 
and  slowing  down— It  is  important  to  note  that  the  criterion 
for  co.,. pari  son  is  the  ratio  of  either  to  the  sum  of  the  two. 
That  is,  the  total  probability  of  one  or  the  other  happening 
(excluding  leakage)  must  be  unity,  For  this  reason,  de* 
creasing  one  has  the  effect  of  increasing  the  other. 

te  beryllium  moderator  is  the  predominant  con- 
tributor to  the  slowing  down  term  by  virtue  of  its  high 
value  of  ^  •  tr^ny   decrease  in  beryllium  thus  tendc  to  en- 
hance the  absorption  in  uranium  during  slowing  down.  Thus 
the  ratio  of  fast  to  slow  flux  varies  inversely  with  the 
berylliiim  concentration  as  shown  in  Figure  2  where  $*ast/  0th 
is  plotted  as  a  function  of  per  cent  beryllium  for  reactors 
of  various  compositions. 

G  •  F.ffect  of  fuel  and  moderator  concentrations 
On  Tower  Jist. ibution 

As  discussed  above,  either  more  fuel  or  less 
beryllium  increases  the  ratio  of  fast  to  slow  flux  and  thus 
raises  the  ability  of  th«  reflector  to.  return  slowed  neutrons 
to  the  core  and  assist  in  levelling  the  power  distribution. 
Fig.  3  shows  this  effect  for  various  beryllium  concentrations 
in  an  otherwise  identical  reactor.  As  the  beryllium  decreases, 
the  power  near  the  core-reflector  boundary  rises. 
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,e  uniformity  of     .r  distribution,  measured  by 
the  ratio  of  ma     i  to  aver    power,  d<         the  point 
at  which  P_,  _  bcourSo        y  be  either  at  the  center,  the 
edge,  or,  ideally,  both  sluiu     ously.  That  is,       3hould 

to       the  center  power  5r  order  to  increase 
P  Le  withoi        neously  raising  Pffir,y» 

As  an  indication  of  this,  the  ratio  of  P  r    /P    in    .  3 
3  to  1.1      eryllium  decreases  from  30  to  20  per 
cent),  then  riser,  to  1.      'beryllium  is  further  decreased. 
In  the  latter  case  the  ed£f  power  is  13 

Fig.  lir   sho     similar  effect  011  the  power  dis- 
tribution for  various  fuel  concentrations.  ',      »  -!  1  and  4b 
indicate  flux  varJ     is  for  various  fuel  concentration!  „ 
Fig.  4b  sho        with  1  ore  fuel  the  flux 

ults  in        f leant  number  of  reu  returned 

to  the  core  before  bein.     pletely  .  These  aid 

Hi]       lower  distribution. 

Al        ;.  hb   (  .    fuel)  the  large  ratio  of 
st  to  thermal  flux  results  in  some  16>  of  fissions 
takin     ce  in  the  two  higher  en  1     roups,  as   1  trasted 
for  the  reactor  of  Pig.  4a  (O.j       >. 

' C:  *-  Or    C OP'Ki^J . t  \ 

if  feet  of  varying  of  the  three 

reflector  components  \-:?r   inv        d.  Ie  gener  1,  for  the  re- 
flector compositions  studied,  beryllium  makes     largest 
contribution  to  t'-.    Lowing  do     cm  and  the  smallest  to 
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the  absorption.  Therefore,  decreasing  the  reflector's 
beryllium  concentration  effectively  increases  the  absorp- 
tion tendencies  while  decreasing  the  moderation.  Reflectors 
with  less  beryllium  thus  return  fewer  slowed  neutrons  to 
the  core. 

Fig.  5»  which  shows  the  power  distributions  for 
three  reactors  with  identical  cores  but  with  different 
proportions  of  beryllium  in  the  reflector,  demonstrates 
this  effect.  The  rise  in  power  occasioned  by  neutrons 
returned  from  the  reflector  is  greatest  for  the  reactor 
having  the  most  beryllium  in  the  reflector,  least  for  the 
reactor  with  least  reflector  berylliur  . 

e.  reflector  Thickness 

Figs.  6,  7,  and  3  show  the  spatial  distribution 
of  flux  for  reactors  with  reflector  thicknesses  of  10,  20, 
and  50  cm  respectively.  Since  fast  neutrons  which  leak  into 
the  reflector  are  subject  to  a  greater  chance  of  slowing 
down  in  a  thick  reflector,  the  thermal  flux  rises  more 
in  such  a  reflector.  In  turn,  this  causes  a  power  rise  near 
the  core-reflector  interface,  since  more  thermal  neutrons 
are  returned  to  the  core. 

This  provides  a  means  of  increasing  the  thermal 
flux  in  the  outer  core  region  where  the  power  normally  is 
less  than  at  the  reactor  center.  By  a  proper  choice  of 
reflector  this  increased  thermal  flux  in  the  outer  core 
region  can  result  in  a  power  level  as  high  as  the  power 
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at  the  reactor  center.  This  represents  the  closest  possible 

approach  to  a  uniform  power  distribution  for  a  given  homo- 
geneous reactor,  since  the  power  distribution  sags  between 
these  two  maximums. 

Fig.  9  shows  the  effect  of  varying  the  reflector 
thickness  on  the  power  distribution  for  a  given  core.  A  10 
cm  reflector  gives  a  power  peak  at  the  reactor  center. 
Somewhere  between  10  and  20  cm  reflector  thickness  the 
optimum  power  distribution  for  this  reactor  occars.  For  20 
and  50  cm  reflector  thicknesses  a  power  peak  occurs  at 
the  surface  of  the  spherical  core.  With  increasing  reflec- 
tor thickness  a  limiting  value  of  Pmax/^avg  ^s 
approached  asymptotically.  A  plot  of  the  neutron  flux  in 
Fig.  3  for  the  50  cm  reflector  shows  that  the  flux  in  all 
energy  groups  approaches  zero  asymptotically  to  the  axis. 
Since  the  current  density  J  =  -Dd(NV)/dr  (as  given  in 
Chapter  I),  a  nearly  zero  slope  to  the  flux  curve  means 
that  there  are  very  few  neutrons  escaping  from  the  outer 
reflector  surface.  This  Indicates  that  a  further  increase 
in  reflector  thickness  would  have  little  effect  on  the 
power  distribution,  so  that  this  is  essentially  an 
infinite  reflector. 

The  positive  slope  (at  the  core-reflector  inter- 
face) to  the  thermal  flux  curve  means  that  J  is  negative, 
indicating  a  net  flow  of  neutrons  from  the  reflector  into 
the  core.  It  is  noted  that  with  reflectors  of  10  and  20 
cm  thickness  there  is  a  net  leakage  of  group  III  neutrons 
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into  the  reflector.  However,  with  a  50  cm  thick  reflector 
there  is  a  positive  slope  to  the  group  III  flux  at  the 
interface,  indicating  a  net  return  of  group  III  neutrons 
from  the  reflector  to  the  core.  This  tends  further  to  in- 
crease the  power  rise  in  the  outer  core  regions. 

The  significance  of  fissions  in  group  III  can  be 
further  investigated.  Figs.  4a, 4b,  6,  7,  and  8  show  flux 
plots  for  different  reactors.  Histograms  in  these  figures 
show  that  approximately  60%   of  the  total  fissions  take 
place  in  energy  group  III.  Thus,  one  would  expect  that  the 
spatial  power  distribution  would  follow  the  core  flux  plot 
of  group  III  quite  closely.  An  examination  of  the  figures 
mentioned  shows  that  the  power  distribution  does  follow 
the  group  III  flux  very  closely  except  in  the  outer  core 
regions  where  the  thermal  flux  begins  to  rise.  At  this 
point  the  power  curve  also  exhibits  an  upward  trend.  This 
reflects  the  fact  that  the  thermal  neutron  fission  cross 
section  is  approximately  three  and  a  half  times  that  of 
the  average  group  III  cross  section,  so  that  a  small  inc- 
rease in  thermal  flux  has  great  significance. 

Another  trend  which  oan  be  noted  from  the  figures 
cited  is  that  the  group  III  flux  distribution  sags  con- 
tinually as  distance  from  the  center  increases.  Thus 
power  which  depends  primarily  on  group  III  will  sag  also, 
even  though  the  thermal  flux  remains  relatively  flat.  If 
more  moderator  is  added  so  that  the  reactor  approaches 
more  closely  a  thermal  reactor,  the  sag  in  the  thermal 
flux  in  the  outer  regions  will  not  be  compensated  for 
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by  a  return  of  reflector  moderated  neutrons,  since,  with 
relatively  few  fast  neutrons  present  there  will  be  little 
fast  leakage  input  to  the  reflector.  Fig.  4a  illustrates 
this  trend. 

f«  Effect  of  Varying  Na/Fe  Ratio 

In  Fig.  11  is  plotted  power  distributions  for 
three  reactors  identical  except  for  the  sodium  and  stainless 
steel  content  of  the  core.  The  curves  are  very  similar 
except  for  the  outer  core  region. 

Because  Fe  has  a  higher  absorption  cross  section 
it  is  to  be  expected  that  increasing  the  core  Fe  content 
will  tend  to  raise  the  fast  flux  somewhat,  thereby 
increasing  fast  leakage  and  the  number  of  slowed  neutrons 
returned  from  the  reflector.  This  point  is  supported  by 
the  plot  of  run  (11)  which  has  the  highest  value  of  power 
at  the  corc-r  f lector  interface,  and  the  highest  Fe  con- 
tent. 

An  additional  effect  is  the  variation  in  diffusion 
coefficient  with  varying  Na/Fe  ratio  in  the  core.  Since  Fe 
has  a  higher  scattering  cross  section  (by  about  a  factor 
of  three)  than  sodium,  and  a  higher  atomic  density  (by  a 
similar  factor)  a  greater  Fe  content  tends  to  decrease  D 
and  thereby  decrease  leakage.  It  is  thought  that  this 
partially  explains  the  slight  difference  between  the  plots 
of  runs  (1)  and  (10)  as  shown  in  Fig.  11.  Run  (1),  with 
less  iron  and  more  leakage  experiences  a  slightly  greater 
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rise  near  the  core-reflector  boundary. 

Since  the  two  effects  mentioned  are  in  opposi- 
tion the  above  discussion  is  rather  hypothetical.  In 
General,  varying  the  Ila/Fe  ratio  has  little  effect  on  the 
overall  power  distribution. 

6«  Optimum  Power  Distribution 

In  Fig.  10  are  plotted  two  optimum  power  distri- 
butions obtained.  Shapes  of  the  curves,  even  though  quite 
different  core  compositions  were  considere  ,  are  practi- 
cally identical.  .This  tends  to  indicate  that  a  limiting 
value  of  power  uniformity  has  been  met.  In  other  words,  the 
power  sag  observed  evidentally  cannot  be  further  lessened 
for  the  reactor  type  and  size  considered. 

The  two  reactors  considered  in  Fig.  10  have 
critical  radii  of  42  cm  and  60  cm.  Since  reactor  size 
affects  the  proportion  of  the  core  which  may  be  affected 
by  reflector  moderated  neutrons,  it  is  thought  that  size 
will  play  an  important  part  in  determining  uniformity  of 
power.  A  reactor  with  a  smaller  core  will  have  a  larger 
proportion  of  volume  which  experiences  a  power  rise  due  to 
reflector  moderated  neutrons.  However  it  can  be  seen  in 
Fig.  10  that  decreasing  the  size  from  60  to  42  cm  has 
essentially  no  effect  on  power  distribution. 

The  two  reactors  considered  i    I  .  10  have 
different  material  compositions  and  different  critical 
radii.  They  are  representative  of  the  class  of  reactors 
investigated.  It  can  be  concluded,  therefore,  that  in 
the  range  of  size  and  composition  investigated  the  optimum 
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power  carve  has  a  characteristic  shape  which  is  that 
shown  in  Fig.  10.  The  ratio  of  PmaxAavg  for  this  curve  is 
approximately  1.16. 

h .  Comparison  of  Four  and  Six  G-roup  ..ethods 
Fig.  13  gives  a  comparison  of  calculations  by 
four  group  and  six  group  methods  for  the  same  reactor. 
In  general  the  results  are  in  good  agreement,  since  the 
critical  radii  differ  by  less  than  four  per  cent.  The  ratio 
of  fast  to  thermal  flux  was  3.4  for  four  groups  and  4.7 
for  six  groups.  These  values  are  of  the  same  order  of  magni- 
tude. 

The  plot  of  spatial  power  distribution  in  the 
core  shown  in  Fig.  13  gives  different  values  of  Pmax/pavg 
for  the  two  methods,  although  in  general  the  curves  have 
similar  characteristics*  ""otter  agreement  in  power  dis- 
tribution between  methods  seemed  to  exist  for  reactors 
having  other  compositions,  although  other  runs  o'n  identi- 
cal reactors  were  not  made  due  to  time  limitations.  Since 
both  methods  use  the  age  diffusion  equation  (with  minor 
differences  which  are  covered  in  the  discussion  of  results) 
the  results  should  be  very  similar.  Experimental  results 
for  Na,  Be,  Fe  -intermediate  reactors  were  not  available 
for  a  check  on  the  accuracy  of  the  results. 
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discussio;:  of  results 

e  principal  results  of  this  study  have  been  the 
determi  atlon  of  critical  radius,  neutron  flux  and  fission 
power  distribution,  and  the  effect  upon  these  quantities  of 
varying  the  material  compositions  of  the  core  and  reflector 
as  well  as  varying  the  reflector  thickness. 

e  method  of  attack  involved  modifying  the  age 
diffusion  equation  so  as  to  permit  solutions  to  be  ob- 
tained by  use  of  Whirlwind  I  digital  computer.  One  method 
was  developed  in  Chapter  1  which  used  a  finite  difference 
approximation  to  the  age  diffusion  equation  with  neutrons 
divided  into  four  energy  groups.  A  second  method  was  devel- 
oped in  Chapter  II  which  used  a  method  of  numerical  inte- 
gration with  neutrons  divided  into  six  energy  groups. 

...e  results  of  calculations  based  upon  Identical 
reactors  using  the  two  different  methods  were  in  substantial 
agreement.  !  owever,  minor  differences  in  spatial  power  dis- 
tributions and  critical  radius  were  noted  for  solutions  by 
the  two  methods. 

The  explanation  for  differences  in  the  radius  and 
power  distribution  given  by  the  two  methods  is  thought  to 
be  the  result  of  differences  in  calculated  average  cross 
sections  for  the  various  energy  groups.  The  energy  groups 
are  so  wide  that  it  is  necessary  to  obtain  an  average  cross 
section  and  diffusion  coefficient  for  each  energy  group.  1 
process  of  averaging  could  result  in  incompatibilities 
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between  aver-    aterial  properties  calculated  for  the 

two  methods,  particularly  in  regions  where  the  cross  sec- 
tion changes  rapidly  with  energy,  as  is  the  case  for  the 
scattering  resonance  of  sodium  near  2700  ev. 

Another  major  difference  "between  the  methods 
involved  the  determination  of  the  elastic  scattering 
input  from  energy  group  i  to  the  next  lower  energy 
group  i+1.  In  the  four-group  method  a  factor  f  was  chosen 
such  that  Qiaw=  f  q.    .  The  solution  was  not  sensitive  to 
the  value  chosen  for  f  as  shown  by  successive  runs  using 
improved  values  of  f .  Chapter  I  (page  4)  discusses  the 
manner  in  which  f  can  be  improved  from  run  to  run  on  a 
trial  and  error  basis.  The  use  of  the  f  factor  gives  good 
results  for  wide  ener  y  groups.  The  disadvantage  of 
having  to  estimate  f  can  be  overcome  by  considering  a 
greater  number  of  energy  groups  where  the  variation  of 
the  slowing  down  density,  q,  over  the  group  can  be  consi- 
dered to  be  linear.  In  this  case  sufficient  accuracy 

,,  in   out, 
could  be  obtained  by  the  assumption  that  q«    -  e(q^  +q±      ). 

avg     x 

In  the  six  group  method,  however,  the  assumption 
is  that  the  elastic  scattering  out  of  an  energy  group  is 


equal  to  the  average  value  of  the  slowing  down  density 

•avg" 


in  an  energy  group,  q*        =  q.   .  This  is  eouivalent  to 

xa ".     1 


assuming  that  the  f  factor  ecuais  unity.  These  two 
different  assumptions  concerning  the  elastic  scattering 
out  of  an  energy  group  probably  account  principally  for 
any  differences  between  results  obtained  by  the  two 
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metnods.  This  difference  can  be  illustrated  by  the  fol- 
lowing' figure j 

Core  Slo ;  Do.m  >.>vj  ■-■ity 

vs 
Lethargy 

^qi  avg. 

OUT,  _ZL-avg 


Solid  line  represents 


o 


4.606   10.81 

Lethargy,  u 


"TO'^  (thermal) 


The  solid  line  represents  the  true  distribution  of  slowing 
down  density.  The  slowing  down  density,  q,  increases  for  low 
lethargy  values  corresponding  to  the  birth  energies  (10  to 
0.1  Hev)  of  fission  neutrons.  For  lethargies  greater  than 
the  birth  lethargy,  q  decreases  due  to  absorption  and 
leakage.  For  the  relatively  wide  lethargy  groups  shown  it 
can  be  seen  that  the  assumption  used  with  the  six  group  method 
that  qfVg  =  o°u  is  not  nearly  as  true  as  the  representation 
of  the  true  continuous  slowing  down  curve  by  the  assumption 
that  qfVg  =  f  qjut  with  a  value  of  f  chosen  so  as  to 
approximate  the  true  curve  as  closely  as  possible. 

The  disadvantages  of  using  a  relatively  small 
number  of  energy  groups  are  as  follows:  (1)  the  true 
elastic  scattering  out  of  an  energy  group  is  difficult  to 
estimate  correctly;  (2)  determining  average  neutron  cross 
sections  for  a  particular  lethargy  group  introduces  errors, 
particularly  when  cross  sections  vary  rapidly  with  energ;  . 
Considering  a  larger  number  of  groups  greatly  reduces 
these  problems. 

The  following  advantages  result  from  using  a 
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small  number  of  energy  groups:  (1)  it  is  possible  to 

neglect  elastic  scattering  across  groups  (for-  the  materials 
considered);  (2)  the  energy  spectrum  of  fission  neutrons 
need  not  be  known  with  great  accuracy  if  group  widths  are 
such  that  all  fission  neutrons  can  be  considered  to  be  born 
in  one  or  two  energy  groups, 

ie  four  group  program  is  arranged  to  include  the 
effects  of  inelastic  scattering.  However,  in  this  inves- 
tigation inelastic  scattering  is  neglected,  since  it  is 
considered  to  be  of  only  minor  importance  for  group  widths 
and  materials  considered. 

The  plots  of  neutron  flux  and  power  distribution 
gave  consistent  results  indicating  qualitative  trends  as 
parameters  were  varied.  The  parameters  consider-"'  were: 
(1)  fuel  concentration;  (2)  reflector  composition;  (3)  Na/3e 
ratio;  (4)  Na/Fe  ratio;  and  (5)  reflector  thickness.  An 
effort  was  made  to  choose  reactors  which  would  have  the 
maximum  power  occurring  equally  at  the  reactor  edge  and 
center  to  get  the  smallest  value  of  Pmax/^avg  (ideally 
equal  to  unity).  The  best  value  that  could  be  obtained  was 
1.16  for  run  (8).  Fig.  11  shows  plots  for  two  optimum 
power  distributions.  Comparison  of  these  two  curves  with 
other  power  distributions  near  the  optimum  seems  to  indi- 
cate that  the  general  shape  of  the  curve  cannot  be  changed 

very  much  for  the  range  of  parameters  investigated.  However, 
these  results  would  not  apply  to  reactors  having  a  much  larger  or 
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smaller  critical  radius  since  the  shape  of  the  optimum 
power  curve  would  very  likely  be  considerably  changed  with  marked 
changes  in  reactor  size. 

A  ratio  of  1.16  for  Fff.ax/ravg  wac  the  best,  value 
obtained  for  the  reactor  size  and  type  discussed. 

The  two  programs  which  were  developed  for  use  with 
the  Whirlwind  I  digital  computer  operated  very  satisfactorily, 
and  it  is  hoped  that  other  investigators  will  find  these 
workable  programs  to  be  helpful.  Both  programs  were  character- 
ized by  rapid  convergence  to  the  answer,  averaging  only 
about  four  or  five  iterations  before  the  convergence  criteria 
were  met. 

For  further  work  on  the  subject  of  uniform  power 
distribution,  it  is  L.t_  ested  that  multi-region  reactors  be 
investigated.  The  Whirlwind  programs  developed  in  Chapters  I 
and  II  can  be  modified  so  as  to  apply  to  reactors  having 
nev^ral  regions  of  uniform  material  properties.  The  programs 
would  lengthen  considerably,  however,  and  it  probably  would 
be  necessary  to  use  auxiliary  Whirlwind  storage. 

It  would  be  of  interest  to  compare  the  results  of 
this  theoretical  investigation  with  experimental  results  in 
order  to  determine  how  closely  a  finite  difference  or 
numerical  integration  approximation  to  the  age  diffusion 
equation  will  represent  reality. 

An  extension  of  this  study  to  investigate  reactors 
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of  a  different  range  in  size  using  other  moderators 
and  coolants  would  supplement  the  results  given  here. 
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SUMKARY. 
It  was  the  objective  of  this  thesis  to  study 
the  practicability  of  achieving  a  uniform  power  distri- 
bution in  a  homogeneous  nuclear  reactor,  through  return 
of  moderated  neutrons  from  a  reflector  to  the  core  to 
compensate  for  leakage*.  ;:ulti-group  methods  of  solution 
of  neutron  diffusion  equations  were  developed.  While  it 

was  impossible  to  compare  the  theoretical  results  obtained 

» 

with  experimantal  data,  because  of  lack  of  such  data,  the 
results  appear  reasonable  in  that  they  demonstrate  trends 
which  are  in  agreement  with  physical  considerations. 

Studies  were  based  on  a  single  basic  reactor 
chosen  because  of  suitability  to  high  output  mobile 
applications. 

It  was  found,  for  the  reactor  type  and  size  studied 
tnat,  wnlle  complete  uniformity  of  power  could  not  be 
achieved,  nevertheless  it  was  possible  to  approach 
uniformity  quite  closely.  The  best  value  obtained  for 

pmax/pavg  was  1'16' 
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APPENDIX  A 

For  a  typical  power  reactor,  Olasrtone  gives  an 

estimated  volumetric  power  output  for  the  core  region  of 

3 
600  watts/cm.  •   It  will  be  of  interest  to  determine  the  ap- 
proximate power  potential  of  the  reactor  of  Hun  8,  vhioh 
gave  optimum  power  distribution  with  the  maximum  power 
points  equal  at  the  core  edge  and  center.   The  critical  core 
radius  was  60.00  cm.  with  sodium  comprising  70  per  cent  of 
the  core  volume. 

Because  of  time  limitations,  detailed  heat  transfer 
calculations  are  not  made.   However,  an  estimate  of  the 
power  potential  of  the  reactor  can  be  made  based  upon  the 
following  calculations; 

Core  radius  -  60.00  cm.  «  1.97  ft. 

Core  volume  -  9*05  x  10-5  cnP 

Total  power  output  at  an  assumed  600  watts/cm^ 
Z  5^3  KV 

It  will  be  assumed  that  the  sodium  coolant  speed  is 

5  ft. /sec.   For  a  mean  temperature  of  800°F,  the  sodium  has 

2 

the  following  properties  : 

P  Z  52  lb/ft3 

c  =0.31  Btu/lb°F 

1 — 3.  Olaa stone,  "Principles  of  Nuclear  Reactor  Engineering  ," 
D.  Van  No3trand  Co.,  Inc.,  Mew  York,  1956,  page  636. 

2 — "Reactor  Handbook:   Engineering , "  U.  S.  Atomic  Energy 
Commission,  August,  1955 >    (Oeneva),  page  266. 
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For  power  applications,  it  is  hypothesized  that  a 
cylindrical  reactor  of  1.97  ft.  radius  can  be  chosen  which 
will  have  a  similar  power  distribution  and  volume  to  the 
spherical  reactor.  For  5^3  MW  power  output,  the  mean  tem- 
perature rise,  AT,  of  the  coolant  flowing  axially  through 

70%   of  the  reactor  cross  section  is: 

£  Btu 
T  _  5^-3  MW  x   3.M3  x  10b  Hr~MW 

*   feT  *  3600^  x  52  ^J  x  o.3l  &f   *  0.7  *   (1.97)2ft2 


ft 


=  5^6°F 


Thus  a  power  output  of  5^3  MW  can  be  obtained  from 
this  reactor  if  sufficient  heat  transfer  surface  exists  to 
give  a  mean  temperature  rise  of  5^6  F  with  a  sodium  coolant 
velocity  of  5  ft/sec.   With  an  inlet  temperature  of  600  F 
and  a  mean  outlet  temperature  of  1150°F,  the  maximum  tem- 
perature  of  a  thin  fuel  element  should  not  be  excessive. 

The  ability  to  transfer  the  required  heat  from  the 

fuel  elements  to  the  coolant  is  the  primary  limiting  factor 

» 

which  must  be  studied  in  determining  the  maximum  power  output. 

It  is  of  interest  to  cocnare  the  output  of  this  re- 
actor with  the  power  requirements  of  a  B-32  type  Jet  botober. 
At  563  miles  per  hour,  eight  Jets  of  10,000  lbs.  thrust  each 
develop  a  total  of  120,000  thrust  horsepower.   If  an  overall 
efficiency  of  20$  is  assumed,  this  is  equivalent  to  a  heat 
source  developing  ^50  M*^°   This  can  be  compared  to  the 
roughly  estimated  reactor  power  output  of  5^3  MW. 
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appendix  b 
table  i 
ccroscopig  gross  sections* 

A«  Scattering  Cross  Sections,  cTs   (bams) 


Energy 

Lethargy,    u= 

E      (Mev) 

ln(10/E) 
0 

Fe 

Be 

U-235 

6.55 

0 

10 

2.4 

3.0 

1.7 

1.25 

.62 

0." 

2.4 

3.6 

1. 

6.45 

1.5 

3.16 

1.151 

2.4 

3.3 

2.1 

6.45 

1.5 

1.778 

1.727 

2. 

2.9 

2.5 

5.24 

1.7 

1 

2.303 

3.5 

2.8 

2.9 

5.08 

4.0 

0.562 

. 

4.0 

3.0 

3.5 

5.78 

3.5 

0.316 

3.454 

.2 

3.0 

4.0 

7.24 

4.0 

0.1778 

4.030 

4.0 

3.5 

4.9 

.7 

3.5 

0.1 

,4.606 

3.2 

4.C 

5.7 

9.76 

3.5 

60.01    x   10"     5-12; 

4.7 

3.6 

10.2 

3.6 

36.21 

5.640 

4.1 

6.4 

5.9 

10.7 

3.7 

20.90 

6.157 

4.2 

2.4 

10.9 

3.3 

12.70 

6.674 

4.9 

4.7 

5.9 

10.9 

3.5 

7.59       ' 

7.191 

7.7 

8.6 

11.1 

3.8 

4.55 

7. 70S 

17.5 

6.2 

6.0 

10.7 

3.5 

2.69 

c.225 

200.0 

6.8 

10.9 

3.3 

1.65 

8.742 

10.0 

7.3 

6.0 

11.7 

3.8 

0.955     ' 

9.259 

4.5 

•  3 

14.6 

3.5 

C .  569 

..776 

3.2 

10.5 

6.0 

18.2 

3.3 

0.346 

10.293 

3.1 

11.0 

20.6 

3.6 

0 .  200      ' 

'     10 .  810 

"b  12.142 

3.1 

11.4 

6.0 

19.2 

3.8 

55.0  x  10 

3.1 

11.4 

6.0 

17.6 

3.8 

14.45 

'      13.474 

3.1 

11.4 

6.0 

19.0 

3.3 

4.07 

'      14.806 

3.2 

11.4 

6.0 

19.0 

3.8 

1.00 

'      16.133 

3.2 

11.4 

6.0 

3.6 

3.8 

0.269 

'      17.470 

3.2 

11.4 

5.9 

2.5 

3.3 

0.070     ' 

1      13.800 

3.2 

11.4 

5.9 

10.0 

3.8 

lieutron  Cross  Sections,  BNL-325,  'U.S.  Atomic  Energy 
Commission,  Aug.  1955  (Oeneva  Supplement) 
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MICROBCOPIG  CROSS  SECTIONS  (Continued) 


3.  Absorption  and  Fission  Cross  Sections,  o^a   and  &%      (barns) 


Energy  Lethargy,  u« 

E   (j>ev)    ln(10/E) 


U-235 


10 

0 

5.62 

0.576 

3.16 

1.151 

1.778 

1.727 

1.000 

2.303 

0.562 

2.879 

0.316 

3.454 

0.1778 

4.030 

0.1 

4.606 

60.01 x 

10- 

3   5.123 

36.21 

H 

5.640 

20.90 

II 

6.157 

12.70 

II 

6.674 

7.59 

II 

7.191 

4.55 

II 

7.706 

2.69 

11 

8.225 

1.65 

II 

8.742 

0.955 

N 

9.259 

0.569 

II 

9.776 

0.346 

II 

10.293 

0.200 

H 

10 .  310 

55.0  x 10" 

12.142 

14.45 

II 

13.474 

4.07 

M 

14.806 

1.00 

11 

16.133 

0.269 

II 

17.470 

.0.070 

II 

18.800 

Ka 


0 


C 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


i"  e 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.1   0.1 

0.15  0.1 


^e 

orf 

era 

0 

1.29 

1.55 

1.29 

1.55 

0 

1.29 

1.55 

1.30 

1.56 

0 

1.27 

1.52 

1.20 

1.44 

0 

1.30 

1.56 

1.50 

1.80 

0 

1.70 

2.04 

2.00 

2.40 

0 

2.35 

2.82 

2.80 

3.36 

0 

3.4 

4.08 

4.1 

4.92 

0 

5.0 

6.0 

6.0 

7.2 

0 

7.4 

3.9 

9.1 

10.9 

0 

11.5 

13.3 

14.5 

17.  A 

0 

19.0 

22.3 

52.0 

62.4 

0 

55. 

66.0 

30.0 

36.0 

c 

72.0 

36.4 

0 

150. c  : 

L&7.5 

.  .;06 

320    375 

0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.0001 
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TABLE  II 
AVER,     ICROSCOFIC  GROSS  SECTIONS* 
PER  ENERGY  GROUP 

A.  Four  Energy  Groups  (cross  sections  in  bams) 


cr« 


tr 


<y 


a 


o' 


J/ 


I 


Group   I 

(Lethargy 

range  0 

to   4.; 

..a 

3e 

Fe 

0 

U-5 

3.10 

3.175 

3.20 

2.76 

6.64 

3.01 
2.94 
3.16 
2.64 
6.64 

0 
0 
0 
0 
1.597 

0 
0 
0 
0 
1.331 

2.46 

0.0845 
0.209 

0.0359 

0.120 

0.0085 

Group  II 

(Lethargy 

■   range  } 

r.606   to  10.81) 

Na 

Be 
Fe 
0 
u-5 

22.22 

5.94 

7.10 

3.76 

12.90 

21.6 
5.50 
7.01 

3.60 
12.90 

0 
0 
0 
0 
7.85 

0 
0 
0 
0 

c.54 

2.46 

0.0845 

0.209 

0.0359 

0.120 

0.0085 

Group   III    (Lethargy  rarifte 

10.81  to 

18.80) 

::a 

Be 
Fe 
0 
u-5 

3.16 
5.97 

11.4 
3.8 

13.0 

3.07 
5.54 

11.26 
3.64 

13.0 

0.067 

0 
0.083 

0 
106.2 

0 
0 
0 
0 

.1 

2.46 

0.0345 
0.209 
0.0359 
0.120 

Group  T 

(Lethargy 

'     .      ,   E  =  0.070 

ev,    1000  F    thermal 

ila 
2e 
Fe 
U-5 

3.2 

5.9 

11.4 
10.0 

3.11 
5.46 
11.26 
10.0 

C.30 

o.ooe 

0.60 

375. 

0 
i      0 

0 
320. 

2.46 

0.0845 
0.209 
0.0359 
0.0085 

* 


Averaged  with  respect  to  lethargy. 
U-5  stands  for  U235 
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CROSS  SECTIONS  ( Continued] 
B.  Jix  Energy  Groups 


6 (fast):  0.234 

mev 

to 

10  Mev 

5      :    5.5 

kev 

to  0.1 

^34  Mev 

4      :    128 

ev 

to 

5.5  kev 

3      :     3 

ev 

to 

L22  ev 

2      :    .07 

ev 

to 

3   ev 

l(th)   :    .07 

ev 

Na  atomic  den 3 it 

Be 

Fe 

0  (in  U02) 

U-235  (90%, in  UO 

U-23^ 

y  (xlO2^) 

2) 

:  .0254 
•  .12 
.  j84S 

:  ,0494 
.02224 

.00251 

Mo 


.029 

.074 

.0119 

.0/.16 

.00283 

.0026 


45 

.  J9 

.03  59 
.12 

.  .  JS5 

.o.j  04 


Absorption  Cross-iection3  (barns) (includes  fission  absorption) 


Energy  Group 

1 

2 

3 

4 

5 

6 


Na  Be  Fe  0  U-235 

0.3  0.0  u.6  0.0  374.0 

0.1  U.O  0.2  j.O  155. u 

o.j  0.0  0.0  0.0  54. 

0.0  .0  ,.  0.  J  12.4 

0.0  U.J  0.0  u.J  3.1 

.0  0.0  0.0  0.0  1.6 


stic  scattering  uross-Jections  (barns): 


1 

2 

3 

4 
5 
6 


3.3  5.'     11.4  3.9 

3.3  5.Q  11.4  3.8 

3.1  6.0  11.4  3.8 

27.4  6.0     6.4  3.8 

3.6  5.1     4.3  3.5 


2.7      2.2      2.7 
Fission  Cross-iections    (barns): 


1 
2 

3 

4 
5 
6 


11.0 

10.5 

.5 

.3 


320, 

12' 

45.7 

10.4 

2.6 

1.3 


U-23B 

1.6 

0. 

o.5 

0.5 

0.5 

0.5 


8.4 
9.8 
9  .0 

. 
9.0 
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CROSS  SECTIONS  (Continued) 
Since  BNL-325  lists  total  cross  sections  for 
elements  other  than  fissionable  materials  certain 
assumptions  are  necessary  in  order  to  compute  scattering 
and  absorption  cross  sections.  The  following  assumptions 
were  made  with  respect  to  all  elements  other  than 
uranium': 

(1)  the  thermal  absorption  cross  sections  listed 
in  BNL-325  for  2200  m/s  neutrons  was  decreased 
by  a  l/v  relation  to  the  .07  ev  thermal 
energy  considered  here; 

(2)  the  scattering  cross  section  was  assumed  to 
be  the  total  cross  section  at  the  point 
where  the  absorption  cross  section,  scaled 
by  l/v  became  negligible. 

These  assumptions  resulted  in  absorption  in  only  the 
thermal  and  next  higher  groups. 

For  the  uranium  the  fission  cross  section  is  listed 
in  BNL-325*  The  absorption  cross  section  was  synthesized 
from  this  under  the  assumption  that  the  capture  to  fission 
ratio  is  constant.  The  scattering  cross  section  of  u-235, 
and  all  cross  sections  of  U-238,  proved  to  be  unimportant 
for'  this  problem.  The  capture  to  fission  ratio  of  u-235 
was  taken  to  be  0.184. 
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T  A3  LLC  HE 
ICTOR   COKPOSITK 


Volume      Per     Cent 


Core 

Reflector 

Run  No. 

Na 

Pe 

Fe 

"uOp 
0.3 

Na 

Re 

Fe 

1 

60 

30 

9.7 

20 

75 

5 

2 

60 

30 

9.6 

0.4 

20 

75 

5 

3 

60 

30 

9.5 

0.5 

20 

75 

5 

4 

70 

20 

9.6 

0.4 

10 

85 

r 

5 

70 

20 

9.6 

0.4 

20 

75 

5 

6 

70 

20 

9.6 

0.4 

40 

55 

5 

7 

60 

30 

9.7 

0.3 

10 

85 

5 

8 

70 

20 

9.7 

0.3 

10 

85 

5 

9 

80 

10 

9.7 

0.3 

10 

35 

5 

10 

50 

30 

19-7 

0.3 

20 

75 

5 

11 

40 

30 

29.7 

0.3 

20 

75 

12 

75 

15 

9.8 

0.2 

5 

90 

5 

13 

70 

20 

9.6 

0.4 

15 

80 

5 

14 

60 

30 

9.8 

0.2 

20 

75 

5 

15 

60 

30 

9.8 

0.2 

40 

55 

5 

16 

30 

60 

9.8 

0.2 

20 

75 

5 

17 

30 

60 

9.8 

0.2 

40 

55 

5 

18 

90 

0 

9.7 

0.3 

25 

70 

5 

19 

60 

30 

9.7 

0.3 

2. 

■  5       95 

2.5 
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TABLE  Ig 

ID  INPUT  DATA 
FOR  /OUR  ENERGY  GROUPS 


Run  :.o. 

Energy 
Group 

I 

D 

f 

/2a 

0.02825 

^a 

i^f 

1(C)* 

1.867 

0.756 

0.000106 

0.000219 

II 

0.5686 

1.412 

0.07548 

0.000523 

0.00107 

III 

0.9769 

1.44 

0.05243 

0.003781 

0.01444 

T 

0.9868 

1.00 

C.5196 

0.03471 

0.05247 

1(R) 

I 

1.138 

0.68 

0.06154 

0.0 

II 

0.531 

1.08 

0.12218 

0.0 

III 

0.593 

1.05 

0.11539 

0.000692 

T 

0.601 

1.00 

0.10581 

0.00566 

2(C) 

I 

1.671 

0.50 

0.02819 

0.0001412 

0.000291 

II 

0.569 

1.75 

0.07540 

0.000697 

0.00143 

III 

0.983 

1.28 

0.05233 

0.011136 

0.01929 

T 

0.990 

1.00 

0.05186 

CO  429  8 

0.069 69 

2(R) 

I 

1.138 

0.56 

0.06154 

CO 

II 

0.526 

1.36 

0.12236 

CO 

III 

0.593 

1.18 

0.11539 

0.000692 

T 

0.601 

1.00 

0.10581 

0.00460 

3(c) 

I 

1,870 

0.50 

0.02627 

0.0001768 

0.000364 

II 

0.567 

2.50 

0.07548 

0.000871 

0.001790 

III 

0.982 

1.57 

0.05242 

0.013490 

0.02412 

T 

0.986 

1.00 

0.05194 

0.051228 

CO  6760 

3(R) 

I 

1.138 

0.56 

0.06154 

CO 

II 

0.531 

1.15 

0.12218 

0.0 

III 

0.593 

1.05 

0.11539 

0.000692 

T 

0.601 

1.00 

0.10581 

coo  566 

4(C) 

I 

2.215 

0.621 

0.21473 

0.0001412 

0.000292 

II 

0.580 

1.528 

0.065327 

0.000697 

O.OC143 

III 

1.169 

1.660 

0.038105 

0.011306 

0.01929 

T 

1.194 

l.OC 

0.037816 

0.043668 

0.06939 

*(C)  stands  for  reactor  core  property 
properties. 


;s,  (R)  ."or  reflector 
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WHIRLWIND   INPUT   DATA 


(continued) 

Run 

No. 

Energy 
Grout) 

I 

II 

III 

T 

D 

1.039 
0.522 

0.537 
0.544 

f 

0.68 
1.08 

1.05 
1.00 

0.06884 

0.13231 
0.12968 
0.12820 

£. 

»u 

4R 

0 

0 

0.000522 

0.003914 

50 

I 

II 

III 

T 

2.215 
0.580 
1.189 
1.194 

0.44 
2.50 

1.57 
1.00 

0.021473 
0.065327 
0.038105 
O.O37816 

0.0001412 

0.000697 
0.011306 
0.04367 

0.000292 
0.00143 
0.01929 
0.06989 

5R 

I 

II 

III 

T 

1.138 
0.531 
0.593 
0.601 

0.57 
2.50 
2.00 
1.00 

0.06154 
0.12218 
0.11539 
0.10581 

0 

0 
0.000692 
0.00566 

6C 

I 

II 

III 

T 

2.215 
0.580 
1.189 
1.194 

0.5 

2.0 

1.57 
1.00 

0.021473 
0.065327 
0.038105 
0.037816 

0.0001412 

0.000697 
0.011306 
O.O43668 

0.000292 

0.00143 
0.01929 
0.06989 

6r 

I 

II 

III 

T 

1.4004 
0.5445 
0.7501 
0.7581 

0.85 
1.22 

1.02 
1.00 

0.04694 
0.10209 
0.0868 
0.08587 

0 
0 

0.001033 
0.006 

?o" 

I 

II 

III 

T 

1.867 
0.568 

0.9797 
0.9868 

0.621 
1.528 
1.660 
1.00 

0.02825 
0.075^8 
0.05243 
0.05196 

0.00C106 

0.000523 
0.008781 

0.03471 

0.000219 
0.00107 
0.01444 
0.05247 

7R 

I 

II 

III 

T 

1.034 
0.520 

0.535 
0.5^0 

0.68 
1.08 

1.05 
1.00 

0.06944 
0.1331 
0.1331 
0.1293 

0 
0 

0.000491 

0.003773 

8C 

I 
II 

III 

in 

I 

2.209 
0.5796 

1.184 
1.190 

0.621 
1.528 
1.660 
1.00 

0.02158 

0.06535 
0.03813 
0.03784 

0.000J06 
0.0005231 
0.005871 
0.035^ 

0.000219 

0.00107 
0.01444 
0.05247 

8R 

I 

II 

III 

T 

1.039 
0.522 

0.537 
0.544 

0.68 
1.08 

1.05 
1.00 

0.06884 
O.I3231 
0.12968 
0.12820 

0 
0 

0.000522 
0.003914 
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WHIRLWIND  INPUT  DATA 
(continued) 


Run 
No. 

Energy 
Group 

I 
II 

III 
T 

D 
2.700 

0.591 
1.498 
1.500 

f 

0.01436 
0.05527 
0.02386 

0.02373 

^a 

2^  ^  f 

9C 

0.621 
1.528 
1.660 
1.00 

0.0001)6 
0.000523 
0.009121 
0.03631 

0.000219 
0.00107 

0.01444 
0.05247 

9R 

I 

II 

III 

1.034 
0.520 

0.535 
0 .  540 

0.68 
1.08 

1.05 
1.00 

0.06944 

0.1331 
0.1311 
0.1293 

0 

0 
0.000491 
0.003773 

LOC 

I 

II 

III 

T 

1.698 
0.564 

0.779 
0.784 

0.9 
1.2 
1.4 

1.0 

0.02919 
0.07287 
0.05522 
0.05474 

0.00011 
0.000523 

0.00931 
0.03904 

0.000219 

0.001073 
0.01444 

0.0524? 

LOR 

I 

II 

III 

T 

1.138 
0.531 
0.593 
0.601 

1.0 

.95 
1.0 
1.0 

0.06154 
0.12218 

0.11539 
0.10581 

0 

0 
0.000692 
0.00566 
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APPEHDIX  G 
Derivation  of  Equations 

Equation  (3«5<1)  of  Chap.  I  was  stated  without  proof 
and  used  as  a  basis  for  the  boundary  conditions  of  the  four 
group  method.  A  proof  of  equation  (3.5d)  follows: 

Prove  that 

as  a  preliminary  step  it  is  necessary  to  prove  the 
following  relationship 


h 


Proof 


™    t=<  +  (r~)('»«-'») 


7    n   v  ^   -r  /\  •"  ' 


H  1  „      i  f~ 


*    ^  >„ 


wf  / 


* 


xf^'~V+iH    -*yi     j 


kfc^+^+Afc'*'*} 
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APPEIIDIX        (continued) 


using  the  relationship 

h        -  h      -Ah 
hi- 1        *» 

we  have 


C 
N 


n-  i 

Thus,  equation  (a)  is  proved.  !3ext  introduce  the  variable  ^P 
and  use  the  definition  of  v  in  spherical  coordinates 


h  h 


use  equation  (b)  In  (a)  with  ^lf   replacing  0  j 
where  \U   =r  -■  "* '    r Trr~~ •   Then 
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AppenQi:-:  G    (continued) 


/  rPty     ! 


=  — ) y  **   w  -  **-/,  a/, 

^  t    >i  =  /  hay 

i 


h 


V+i 


V-l/ 


\       V*h  -zv«  +  V„-, 


h  (vt  -zvt  +%) 


+ 


<VZ-  *v,  +  vj 


^>~ 


I  At*  G 


with     Ah   -      r     ;     -   h 


+  i  *        the  result  becomes    (continued) 


—    r    [\J/       -  V- 


^^ 


J^k     4  ^_^  _  ^Ji     + 


ZA* 


r 


"/ 


4- 


A  ^  ^/" 


+  i%    + 


c 

%-*%  +  % 

c 


j           j     %           £    K^N-t,     *~  ^N'l) 

2  a>~ 

Ah             %  ah             6 
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■•:i:DIX  i  (continued) 


since  ^  is  to  be  identified  with  ^0  ~  °   '   Thls  is  the 
proof  of  equation  (3»5d)  of  Chap.  I  which  is  concerned 
only  with  the  first  part  of  the  above  expression  relating 
to  evaluation  at  the  point  rN« 


Hoover        '«3t)0 
Uniform  power  distri- 
ion  in  a  c  Led 

homogeneous  reactor. 


3G5 


i 

in  a  reflect 
or  • 


